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We show that the Dyson Brownian Motion exhibits local universality after a very short time 
assuming that local rigidity and level repulsion hold. These conditions are verified, hence 
bulk spectral universality is proven, for a large class of Wigner-like matrices, including 
deformed Wigner ensembles and ensembles with non-stochastic variance matrices whose 
limiting densities differ from the Wigner semicircle law. 
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1. Introduction and motivation 

In his groundbreaking paper Wigner conjectured that the eigenvalue gap distribution of large 
random matrices is universal and that it serves as a ubiquitous model for the local spectral statistics 
of many quantum systems. The Gaussian case was fully understood in the subsequent works of Dyson, 
Gaudin and Mehta; see |44| for a summary. This simplest case can be generalized in two directions. 
For invariant ensembles, the joint density function of the eigenvalues can be explicitly expressed in 
terms of a Vandermonde determinant; a formula that can also be interpreted as the Gibbs measure 
of a gas of one-dimensional particles with a logarithmic interaction. For specific values of the inverse 
temperature /3 = 1,2,4, the correlation functions may be expressed and analyzed using asymptotics 
of orthogonal polynomials [31] and universality was proved under various conditions on the potential 
in HZl mi 03 [50], with many consecutive works following. This method, however, is not applicable 
for other values of /3 even in the Gaussian case, where the correlation functions were described in [54] . 
Universality for general /3-ensembles was first established recently in mm for (3 > 1, with different 
proofs given later in 13 EH that also hold for (3 > 0 albeit with more restrictions on the potential. 

Among the non-invariant ensembles, the most prominent case is the N x N symmetric or hermitian 
Wigner matrix characterized by the independence of the entries (up to the constraint imposed by the 
symmetry class). Beyond the Gaussian case there is no explicit formula for the eigenvalue distribution 
in general, but in the hermitian case (/3 = 2) and for distributions with a Gaussian component, the 
correlation functions can still be expressed using an algebraic identity (Harish-Chandra-Itzykson-Zuber 
integral). A rigorous analysis of this approach yielded universality for hermitian Wigner matrices with 
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a substantial Gaussian component, 133 0. The first proof of hermitian Wigner universality for an 
arbitrary smooth distribution was given in [21], the smoothness condition was later removed in [531 [25 j. 
Lacking the algebraic identity, the symmetric case (/3 = 1) required a completely different approach 
based on the analysis of the Dyson Brownian motion. The basic observation of Dyson m was that the 
eigenvalues of a matrix ensemble, embedded in a simple stochastic flow ( Dyson matrix flow), evolve 
autonomously and satisfy a system of N stochastic differential equations, called the Dyson Brownian 
Motion (DBM). The universal eigenvalue statistics emerge in the bulk spectrum as a consequence of 
the invariant measure of the DBM. Local statistics require to understand only the local equilibration 
mechanism which occurs on a very small time scale that can be bridged by perturbative methods. The 
rigorous theory of this idea was initiated in [23] and developed in a series of papers [251 [50] leading 
to the complete proof of the Wigner-Dyson-Mehta universality conjecture for Wigner matrices in all 
symmetry classes; see m for a summary. More recently two stronger versions of the bulk universality 
have been proved. In contrast to the previous results that required a local averaging, the universality 
of each single gap was shown to be universal in [25], while the universality of correlation functions at 
each fixed energy was obtained in m • These papers heavily relied on a new tool from [28], the concept 
of Holder regularity theory for the parabolic equation with random coefficients given by the DBM. 

In all these works on the spectral universality for Wigner matrices, the global limiting density was 
the semicircle law; in particular it did not change in time under the DBM. The same Gaussian measure 
and its localized versions could be used as equilibrium reference measures for all times. The main idea 
was to artificially speed up the global convergence by considering the local relaxation flow [24] and 
then to prove that the additional local relaxation terms do not substantially modify the local statistics 
thanks to a-priori bounds on the location of the particles. These bounds are called rigidity estimates 
and they directly follow from short scale versions of the Wigner semicircle law that are called local laws. 

The method of local relaxation flow has two main limitations that are related. First, it operates with 
global measures, in particular, quite precise rigidity information is needed for all eigenvalues. This is 
clearly unnecessary (and in some cases hard to obtain); far away eigenvalues should not influence local 
statistics too much. Second, if the initial matrix of the Dyson matrix flow does not obey the semicircle 
law, then the density changes with time following the semicircular flow, related to the complex Burgers 
equation for its Stieltjes transform. The time dependence of the density was originally not incorporated 
in the method of the local relaxation flow. This second limitation was tackled very recently in |4flj . 
where universality for deformed Wigner matrices with large diagonal elements was proved. Using ideas 
from hydrodynamic limits [58] . a global reference measure was constructed as an invariant /3-ensemble 
with a “time” parameter so that its equilibrium density trails the semicircular flow. This equilibrium 
measure was then used as a basis to construct the local relaxation flow. Once the fast convergence to 
the reference measure is established one can infer to the universality of the /3-ensemble mm, or, 
alternatively, one can use the uniqueness of the local Gibbs measure established in [28] to conclude 
universality with a tiny Gaussian component. This result is then easily complemented by a standard 
Green function comparison method to remove the Gaussian component entirely. 

As a technical input for the analysis in {40], the global rigidity for the reference /3-ensemble is 
required, which is not available for the case when the equilibrium density is supported on several 
intervals. In particular, the result of m is limited to the deformed Wigner ensembles with a single 
interval support that excludes the case when the diagonal has a strongly bimodal distribution. 

We remark that bulk universality for special classes of deformed Wigner matrices in the hermitian 
symmetry class has also been proven with different methods. The local sine kernel statistics for the 
sum of a GUE and a diagonal matrix with two eigenvalues ±a of equal multiplicity has been obtained 
with Riemann-Hilbert method mmm- In particular, the density in this model is supported on two 
disjoint intervals if a is sufficiently large. The GUE matrix can be replaced with an arbitrary Hermitian 
matrix if the first four moments of its single entry distribution matches those of the Gaussian [35] . 
A much more general class of deformations of the GUE has been tackled in [52] relying on a version 
of the Harish-Chandra-Itzykson-Zuber integral. Using Green function comparison techniques [29] and 
the local laws from [391 EH Eg], one can replace the GUE with any hermitian Wigner matrix under 
the four moment matching condition. 

Random matrices whose limiting densities are supported on several intervals arise in other promi- 
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nent contexts as well. We call symmetric or hermitian matrix ensembles, H = (hij ), Wigner-like if 
their entries are independent (up to the symmetry constraint). If, in addition, the matrix elements are 
centered, E hij = 0, and the sum of the variances Sij = E|/qj| 2 in each row is constant, say one, i.e., 

N 

Y. Sij = const = 1, VI, (1.1) 

3=1 

then the limiting density is the semicircle law. If either condition is violated, the limiting density 
is generally not the semicircle law and typically it may be supported on several intervals. The case 
of H = W+A, where W is a standard Wigner matrix with i.i.d. centered entries and A is a deterministic 
matrix (representing the nonzero expectations E h^), was considered in [35]136]. where local laws and 
rigidity were established. If condition (O) is dropped, then an even richer class of possible limiting 
densities arise. These were extensively analyzed in EJ®, where all possible density shapes are classified, 
local laws and rigidity are proven. 

In the current paper, we prove bulk universality for all these models. As in the previous papers 
using DBM, the key part is to show universality for matrices with a tiny Gaussian component. 

Beyond these applications, our main result is formulated on a more conceptual level. Dyson argued 
in m that the local equilibrium of the DBM is attained after a very short time irrespective of the 
global density. In fact, the global density equilibrates on a time scale of order one, while the local 
equilibration time is of order 1/N. The local equilibration is solely due to the logarithmic interaction 
in the DBM, while the evolution of the global density is given by the semicircular flow. In this paper 
we fully decouple the effects of these two processes. In the main Theorem 12.11 we prove bulk local 
universality for the DBM assuming that it satisfies rigidity and level repulsion, but only locally. On 
the global scale only a very weak version of rigidity is required, in particular the condition is insensitive 
to outliers or to the behavior at the edges. These assumptions can then easily be verified from local 
laws in each model. 

After completing this manuscript, we learned that similar results were obtained independently 
in [38] . 

Notational conventions: We use the symbol 0( ■) and o( •) for the standard big-0 and little-o 
notation. The notations O, o, <C, 3>, refer to the limit N —> oo. Here a <C b means |a| < AT _ ^|6|, 
for some small £ > 0. We use c and C, C' to denote positive constants that do not depend on N. 
Sometimes we use subscripts or superscripts to distinguish ./V-independent constants, e.g., cq, c±, c' etc. 
Their value may change from line to line. Similarly, we will use £ > 0 for a small, respectively D > 0 
for a large positive exponent, mainly appearing in various rigidity bounds. Their precise values are 
immaterial; at the end of the proof it may be chosen sufficiently small, respectively sufficiently large, 
depending on all other exponents along the proof. Finally, we use double brackets to denote index 
sets, i.e., for m,n 2 £ R, 

[ni,n 2 ] := [ni,n 2 ] HZ, N n := [1, Nj . 

Acknowledgement: The authors thank O. Ajanki and T. Kruger for many valuable discussions at 
the early stage of the project. 


2. Main results 

In this section, we give a detailed description of our model, including all assumptions, and state our 
main results. We start with introducing basic concepts such as the Stieltjes transform, the semicircular 
flow and the Dyson Brownian motion (DBM). 

2.1. Stieltjes transform. Given a probability measure, v, on M, define its Stieltjes transform, m v , by 

m v {z) := [ ; z€C + :={z£C,Imz>0}. (2.1) 

Jr v- z 
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Note that m„ is an analytic function in upper half plane. In the following we usually write z = E + ir), 
E G R, 77 > 0, and we refer to E as an “energy” and to z as the spectral parameter. For given 77 > 0, 
we let P v denote the Poisson kernel defined by 


p ^) : =--srr^> p eR, (2-2) 

7T E Z + ifl 

and we note that f R P v {E)dE = 1 and P Vl+V2 (E) = (P m *P V2 ){E), for all 77 , 771,772 > 0, E G 1, where * 
denotes the convolution on R. We further remark that 

^-flmm v {E + i/ 7 ) = (P„ * v){E ). (2.3) 

7T 

Assuming that v admits a density, which we also denote by u, we can recover v from m u through the 
Stieltjes inversion formula 

iy(E) = — lim lmm v (E + irj) = lim(P„ * v){E ), Eg R. (2-4) 

7T ri\0 V\0 

The Hilbert transform, (T v), of v is defined by as the principal value integral 

{Tv) (E) := f ^, EG R. (2.5) 

JrV-E 


2.2. Semicircular flow. We next introduce the semicircular or classical flow. Let A4(R) denote the 
set of probability measures on R. Then the semicircular flow is the process R + x A4(R) —> A4 (R), 
{t, g) 1 —^ Et[g\ obtained via its Stieltjes transform as follows. For t = 0, set Pa[g\ '■= g. For t > 0, let 
rrit(z) satisfy 


m t {z) 


d g{y) 


s t / 2 y — z — (1 — e ( )mt(z) ’ 


Imm t (z) > 0 , z G C + . (2-6) 


It is straightforward to check |48| that (12.61) has indeed a unique solution such that lim inf r/ x 0 Im m t {E+ 
ir]) < 00 , for any E G R, t > 0. In fact, for t > 0, mt has a continuous extension to C + UR [§] that we 
also denote by m*. Set then 

Pt[g]{E) := —lim Im mt {E + ig), t> 0, Eg R, (2-7) 

7r 7 j\,0 


so that Ft[g\ is defined through its density F t [g\{E), Eg R. In particular, for t > 0, Pt\g] is an 
absolutely continuous measure. (For simplicity we use the same symbol for absolutely continuous 
measures and their densities.) 

Further, it is easy to check that m t {z) converges pointwise to 


m 0 {z) 


f dgp(y) 
/r y-z 


( 2 . 8 ) 


for all z G C + , as t \ 0. It follows that Pt[g] converges weakly to g as t \ 0. Starting from (12.61) 
and (12.71) . one also checks that 

Tt+s[g] =E t o p s [g] = F^jFslg ]], t>8, g G A4(R). 

In fact, using the additive free convolution, the flow t J~t can be endowed with a {w* -continuous) 
semigroup structure Esina ma; see also [155 ], [52 ] for reviews. Yet, we will not pursue this point of view 
in the present paper. 

In the following, we often write gt := Pt[g] with go = g and we call t 1 —>■ g t the semicircular 
flow started at g. Recalling m it is clear that rrq is the Stieltjes transform of g t and we simply 
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write mt = m et . We remark that the standard semicircle law, g sc , is invariant under the semicircular 
flow, i.e., Ft[Qsc\ = Qac > f° r all t > 0 , and that g t = Ft[o\ converges weakly to g sc , as t oo, for 
any g G A4(R). This follows directly from (12.61) and the fact that the Stieltjes transform, m Bac = m sc , 
of g sc satisfies m sc (z) = — ( m sc (z ) + z) _1 , z G C + . 

For IV G N and hxed t > 0, let 7 (t) = ( 7 fc(i)) denote the set of iV-quantiles with respect the 
density gt, where 7 k{t) is the smallest number satisfying 

n k(t) 

/ g t (x) dx = — , g t = jF t [g) , (2.9) 

for all t > 0. It is straightforward to check that 7 fe(t) inside the “bulk”, i.e., where g t is strictly positive, 
is a continuous function of t. This follows from the (weak) continuity of the flow t 1 —> g t . Moreover, 
the points 7 (t) in the bulk approximately satisfy a gradient flow of a classical particle system with a 
logarithmic two-body interaction potential between the particles (see Lemma 14.61 below!. We refer to 
Appendix [A] for a more detailed discussion. 


2.3. Dyson Brownian motion. Fix IVgN and let F^ C M. N denote the set 

rW:={A=(A 1 ,A 2 ,...,A JV )cK JV : A x < A 2 < ... < A w } , (2.10) 


and denote its closure by F^ ■ 

Dyson Brownian motion (DBM) is given by the following stochastic differential equation (SDE) 


dAj(t) 



1 

N 


E 


A j(t) - A i(t) At 



* G N n , 


P> 1 , 


( 2 . 11 ) 


with fixed initial condition X(t = 0) G F^ N \ where /3 > 1 is a fixed parameter with the interpretation 
of inverse temperature, and where are a collection of independent standard Brownian motions 

in some probability space (fi, P). We denote by E the expectation with respect to P. 

It is well known, see Section 4.3.1 of [T|, that (] 2 . 1 1 |) with j3 > 1 has a unique strong solution, A(t), 

for any initial condition A(0) G F^- Further, for any t > 0, we have A(t) G F^ almost surely. 

The equilibrium measure for the DBM is the Gaussian invariant ensemble explicitly given by 

1 N I 1 

p G (A)dA^ M W(A)dA=- M e-^dA, U G := £-A? - - £ log (A, - A, ; ), (2.12) 

^/3,G i= 1 1 <i<j<N 

where dA := 1(A G F^) dAidA 2 • • • dAjv and where Z^l is a normalization. For fixed /3, we denote 
by E g the expectation with respect the measure g G in (12.121) . 

Consider next a sequence of vectors A^ (0) = (A^(0),..., A^(0)) G F^ N \ N G N. Let A ^ N \t) = 
(A [ N \t ),..., A ^\t)) G F^ denote the sequence of vectors such that, for each N gN, A l ' N \t) G F^ 
is the solution to (12.111) with initial condition A <Ar l(0). For simplicity we abbreviate A(f) = A (W l(f), 
respectively A i(t) = i G Njv, in the following. 

Assume that there is a probability measure, g q°, on R such that 



i= 1 


as N —> 00 , i.e., the empirical distribution of the initial data A (Ar l(0) converges weakly to pg°. Then, 
under some mild technical assumptions on A (jV l(0), Proposition 4.3.10 of [T] states that 

1 N 

^ N) ■■= ^ QT, t> 0, (2.13) 

V i=1 

as N —> 00 , where g%° := Btig^] denotes the semicircular flow started from c.f., Subsection 12.21 
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2.4. Main result. In this subsection, we state our main result. We need one more definition: A 
labeling £ is a random variable l : R —> Z, x £(x) such that £(x + 1) — £(x) = 1 and £(x) = £(\x\). 

Theorem 2.1. Let A (t), t > 0, be the solution to the DBM in (12.111) with deterministic initial 
condition A(0). Given any small positive e > 0 and any small 6 G [0,1/20], with e > 28, consider times 
, <2 £ with N~ 1+e < t 2 — t\ < N~ e . Let g be a probability measure on R. Denote by g t = T t [g\ 
the semicircular flow started from g. Choose E* G R such that g tl (E*) > c, for some small c > 0. 
Assume that A (t) and g are such that the following conditions are satisfied. 

(1) At time t±, the density g tl = [f?] is regular in the following sense. There is a constant X > 0, 

independent of N, such that the Stieltjes transform m gt of g tl , i.e., 


m etl (z) 


Pti (y) d y 

y-z 


Z G C + , 


(2.14) 


extends to a continuous function on D s := {z = E + vq G C : Eg [E* — X, E* + X], g > 0}, and 
satisfies 

\m eti (z)\<C, \d:m eti (z)\ < C(N 5 ) n , n = l,2, (2.15) 


uniformly on D-%, for some constant C. Moreover, g tl has finite second moment and satisfies 


gti (-®) — c > Eg [E* — X, E * + X], 


(2.16) 


for some c > 0 . 

(2) The process A (t) is rigid and is related to gt in the sense that there is a small a = er(X) > 0, 
independent of N, such that the following holds. 

(a) Strong rigidity inside I a : There is a time-independent labeling £ such that 7 ^( 7 ) (ti) G [ E * — 
X/4, E* + X/4], for all i G I a = \L — crN, L + crlV], where L G is the largest integer such 
that 7 ^(L)(ti) < E *. Moreover, for any (small) £ > 0 and any (large) D > 0 we have 


(jAj(t) - 7 * (i )(i)| < E , Vt G [ti,t 2 ]^ > 1 - N D , Vi G I a 


(2.17) 


for large enough N > Nq(£, D), where (ji(t)) are N-quantiles with respect to the measure gt; 
see (12.91) . 

(b) Weak rigidity outside I a : For any £ G (0,d) and any (large) D > 0, we have 


- E 

N ^ 


1 


gt(x)dx 


k : \L-k\>aN 


A k{t) — E * 


R\I(t) 


x — E* 


< Jfs^tG [tiM ) > 1-N~ D , (2.18) 


(2.19) 


for large enough N > N 0 (£,D), where I(t) := fu{L-aN){t),le{L+a N ){t)]- 

(3) Level repulsion inside I a : For any i G I a and t G [ti, ? 

p(W)-A i±1 (t)| < |A <N s u 0+1 , u> 0, 

for large enough N. 

(4) Holder continuity of DBM: For any (small) £ > 0 and any (large) D > 0, we have 

- Afc (s) | < N^Vt - s , Vt, s G [ti , f 2 ], t > s^j > 1 — N~ d , Vk G Ntv\7ct , (2.20) 

for large enough N > No(fi, D). 

Then, there are small constants f, %, a > 0, such that the following holds. Fix n G N and let 
O : R n —> R be smooth and compactly suppoHed. Fix any T G [t\ + N 2S /N, tf\. Set p* := grirfLiT')')- 
Then 


E 


O [l(Ng*)(\ io (T)-\ io+j 


.7=1 


= E l 


O 


A i' 0 +j 


.7=1 


0(N~ f), (2.21) 
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for N sufficiently large, for any io,i' 0 £ Nat satisfying \iq — L\ < N x , \i' G — L'\ < N x with any 
L' £ [cxN, (1 — a)JV], and where g# := g S c{^lL',sc) denotes the density of the semicircle law g sc at the 
location of the L'-th N-quantile of g sc . 

Remark 2.2. The formula (12.211) expresses the single gap universality, i.e., that the joint distribution 
of n consecutive gaps coincides with that of a Gaussian invariant ensemble for any fixed n. Single 
gap universality clearly implies the weaker averaged gap universality, where (12.211) is averaged over N b 
consecutive ig’s, for some 0 < b < 1. It is well known that averaged gap universality implies the 
averaged energy universality, i.e., the universality of the local correlation functions around an energy E, 
averaged over E near the reference energy E *; see e.g.. Section 7 of [25] , 

Remark 2.3. The measure g in Theorem 12.II is assumed to be deterministic, but it may depend on N 
in contrast to the measure gjf 0 ^ of (12.131) which is indeed the limiting object as N —> oo. Consequently, 
the semicircular flow g t = J~t [g] will also be TV-dependent. Typically one expects that g t converges 
weakly to gf°, yet the speed of convergence may be very slow and hence not be compatible with 
Assumption (2) of Theorem 12. II In Subsection IA.21 we discuss Assumption (1) in more detail. 

Notice that the initial condition A(0) of the DBM and the initial data g of the semicircular flow 
do not have to be related; this will allow us for an additional freedom in the applications. We only 
require that A(f) is close to the quantiles of g t in a short time interval t £ [U, tf\ and only locally near 
the reference energy E *. We also allow for a possible relabeling £ that can be used to accommodate 
outliers in applications. At first reading the reader may ignore t and consider £(*) = i for simplicity. 

2.5. Random matrix flow and universality. In this subsection, we briefly explain how Theo¬ 
rem \m can be used to prove bulk universality for many random matrix ensembles H. We will follow 
the three-step strategy initiated in a series of works [251 [23121]; see m for a concise summary. 

Step 1 is to prove a local law, which includes rigidity for the eigenvalues and bounds on the 
resolvent matrix elements G(z) = (H — z) _1 down almost to the scale of the eigenvalue spacing, i.e., 
for p = Imz TV -1 . This step is typically model dependent, mainly because the limiting density 
of the eigenvalues varies from model to model. The key tool is the self-consistent equation for the 
Stieltjes transform of the density (and its vector version for the individual matrix elements Gu); its 
solvability and stability properties need to be investigated for each model. 

Step 2 is to prove universality for matrices with a small Gaussian component that can conveniently 
be generated by running a matrix valued Ornstein-Uhlenbeck process. Theorem 12.II is used in this step 
and it replaces the previous argument that relied on a global equilibrium measure and its version with 
relaxation. As advertised in the introduction, Theorem [27T| requires rigidity information only locally, in 
particular it also applies to models where the limiting density is supported on several intervals. Step 2 
is model independent once the input conditions of Theorem 12.11 are verified. 

Finally, Step 3 is a perturbation argument which is also very general. Using the Green function 
comparison strategy [29] and the moment matching (introduced first in [53] in the context of random 
matrices), one can remove the tiny Gaussian component. The main input here is the a priori bound 
on the resolvent matrix elements obtained in Step 1. 

More concretely, consider a random N x N hermitian or symmetric matrix H t = Hf with ma¬ 
trix elements ( hij ). Suppose the matrix elements are time-dependent and they satisfy the Ornstein- 
Uhlenbeck (OU) process 


dB 1 

dhij = —j= - -h^ dt, i,j £ Nat , i<j, (2.22) 

where (R,j : i < j) are independent complex Brownian motions with variance t and (Bu) are inde¬ 
pendent real Brownian with variance t for f) = 2; while for /3 = 1, (Bij : i < j) are independent real 
Brownian motions with variance t and Ba are real Brownian motion with variance 2 1. It is easy to 
check that the solution to (12.221) . H t , with initial condition Hq, satisfies the distributional equality 

H t ~ e -t/2 .f/o + (1 — e _t ) 1/2 U, (2.23) 
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where U is Gaussian, i.e., belongs to the GUE (/3 = 2), respectively to the GOE (j3 = 1), and U is 
independent of Hq. 

The eigenvalues of H t , here denoted by A(f), satisfy |T5] the SDE (12.111) . with /? = 1 or /3 = 2, 
where the initial condition A(0) is given by the eigenvalues of the initial matrix Hq. We will run the 
OU process until time t 2 = o(l). Let g denote the limiting density of Hq. We fix an energy E* in the 
bulk spectrum of Hq, i.e., g(E *) > c > 0; it is easy to see that E* stays in the bulk of H t as well for 
any t <t 2 . The assumptions of Theorem 12.11 can then, via the identification (12.231) , be checked from 
the matrix flow H t in the time slice t £ [ti,t 2 ]. The typical choice is t 2 = N~ e and t\ = t 2 — N~ 1+2S , 
with some small positive exponents e -C S. 

Assumption (2) can be checked from a local law for the random matrix H t . We need such informa¬ 
tion not only for the original matrix Hq, but along the whole OU flow. Typically, however, when the 
local law is proven for some matrix Hq, it also holds for H t , i.e., for Hq with a Gaussian convolution. 
Notice that the strong form of rigidity, an almost optimal bound on A i(t) — ”/i(t) expressed in (12.171) . is 
needed only for eigenvalues near E t in the bulk. Much weaker information is needed for far away eigen¬ 
values; the condition (12.181) involves controlling the density only on the macroscopic scale. In terms of 
the Stieltjes transform, m[ N \z) := 77 Tr ( H t — z) _1 , z £ C + , of the empirical density, Assumption (2) 
follows if the bounds 


\m[ N \z) - m t (z)\ < —— , for 2 = E + ig, g £ [N 1+5 ,cE], |E-E*|<E, (2.24) 

A?7 

I m[ N) (E* + ig) - m t (E* + ir?)| < ^ , for g £ [cE, 1], (2.25) 

hold with high probability, for any t £ [ti,< 2 ]- Indeed. (12.251) directly implies (12.18[l . By a simple 
application of the Helffer-Sjostrand formula ( e.g., following the proof of Lemma 8.1 in [20]), we see 
from (12.2411 that 


#{A j(t) £ [in, £2]} - #{7 j(t) G [E U E 2 ]} 


<CN «, \/E 1 ,E 2 £ 


e -\^ e+ ¥ 


( 2 . 26 ) 


In particular, rigidity between the A(f) and 7 (t) sequences holds on scale N~ 1+ % within [E — |E, E + 
^E], This implies |A,;(f) — 7€(»)(^)I < N~ 1+ ^, for any i £ I a , up to an overall shift in the labeling that is 
encoded in the labeling function t(i). We only need to show that the labeling l(i) is time-independent, 
i.e., that along the whole time interval t £ [ti, ^ 2 ] it is the same element of the 7 (f) sequence that 
stays close to a given element of A(t) within the rigidity precision N~ 1+ ^. We call this property the 
persistent trailing of DBM by the flow of the quantiles. Given (12.261) . it is sufficient to check this for 
one element of the sequence; e.g., that if |Al(U) — 7 «z,)(ii)| < N~ 1+ ^ with some shifted index £(L), 
then |Ai(f) — 7 f(z,)(i)| < N~ 1+ for all t £ [t\,t 2 \. Notice that persistent trailing is a nontrivial feature 
of the DBM since the length of the time interval t 2 — t\ = N~ 1+2S is much bigger than the rigidity 
scale N~ 1+ t. Nevertheless, in Proposition IbTTI in Appendix [B] we show that there is an event So in 
the probability space of the Brownian motions with P(So) > 1 — N~ £/ 2 such that 7^(17 (i) persistently 
trails Az,(f). It is easy to see that the universality in Theorem 12.11 also holds if Assumptions (2)-(3) 
are valid only on the event So- 

Level repulsion estimates of the form of Assumption (3) for random matrix ensembles can be 
obtained using the method of [25]. This approach requires two inputs: strong local rigidity as in (12.1711 
and smoothness of the distribution of the matrix elements of H. The former is already verified by 
Assumption (2), the latter needs a slight extension of [25] to “almost-smooth” distributions, where 
smoothing may be provided by the OU process. Indeed, in Appendix B of m it was shown that H t 
satisfies level repulsion in the form (12.191) . if t = N~ cS with some small constant c > 0 (another merit 
of the proof in m is that it also presents the necessary modifications to cover symmetric matrices as 
well, while [25| was written for hermitian matrices only). So we will choose e = in the definition 
t 2 = N~ £ to guarantee that (12.191) holds for any t £ [ti,t 2 \. Notice that the only reason to run the 
DBM up to a relatively large time t 2 = N~ e is to guarantee that the smoothing effect is substantial to 




























yield level repulsion. If the distribution of Hq were smooth initially, so level repulsion in its original 
form |25| applied, we could have chosen t\ = 0, £2 = N~ 1+2S with some small S > 0. 

Finally, Assumption (4) can easily checked as follows. For any two N x N matrices A = A*, 
B = B *, we have dist{Sp(A), Sp(B)} < \\A— B Hoc, where Sp(A), Sp(B) denote the spectra of A, B 
and where || • ||oo denotes the operator norm. Also recall that the operator norm of U is bounded by 
a constant with overwhelming probability; see, e.g., Exercise 2.1.30 of [1]. Thus, choosing A = H t , 
B = H 0 , we see that Assumption (4) is satisfied provided that ||i?o||oo < N^/ 2 with overwhelming 
probability. This bound can be easily proven for all matrix models we have in mind. 

Having checked the assumptions, the conclusion of Theorem 12.II is that gap universality holds for 
any matrix with a substantial Gaussian component of size £2 ~ N~ e . The rest is a standard moment 
matching and Green function comparison argument that we sketch for completeness. 

Given an initial Wigner-like matrix H for which we eventually wish to prove universality, we choose 
t 2 = N~ e with a sufficiently small e > 0. By moment matching (see, e.g., Lemma 6.5 of [29]), we 
construct another matrix Hq such that the solution H t2 at time £2 of the matrix Ornstein-Uhlenbeck 
process (12.221) with initial condition Hq is close to H in the four moment sense. Choosing T = t? in 
Theorem o we obtain gap universality for Ht which also implies universality of local correlation 
functions at E with a small averaging in the energy parameter E around E *. The local eigenvalue 
statistics of H and Ht coincide by the Green function comparison theorem introduced in [29] . More 
precisely, the method of [29] gives coincidence in the sense of correlation functions while Theorem 1.10 
of |54j extends the Green function comparison method to individual eigenvalues, hence to gaps as well. 
This completes our sketch on how to apply Theorem 12.11 for random matrix models. 

2.6. Strategy of the proof of Theorem 12.11 The first part of the proof is to understand the 
dynamics on a macroscopic scale, i.e., to control the semicircular flow and the induced dynamics on 
the time-dependent quantiles 7 *(£). This analysis is of interest itself and it is deferred to the ApoendixlAl 
since it requires quite different tools than the main part of the proof. The key information (collected 
in Section 14.21) is that the quantiles in the bulk move coherently with a local mean velocity that varies 
in time on the macroscopic scale. Since we concentrate on the vicinity of a fixed energy A* and on a 
small time window, by a simple linear shift we can achieve that the mean velocity is negligible near A*. 

The second step is to localize the problem: we choose an integer K 1 such that 

N S K W < TV, K <N S . (2.27) 

We consider the conditional measure on K = 2/i+1 consecutive internal points x = (A l-k, ■ ■ ■, A l+k), 
labeled by I := [A — K, L + KJ, conditioned on the remaining N — K, external points y = (Aj : \i — L\ > 
K). The index L is chosen so that 'ye(L) is close to A*, where 7 *, is the fc-th iV-quantile of the density 
at t\. In the equilibrium setup this corresponds to the local Gibbs measure with boundary conditions 
given by y (this idea was first introduced in [13] in the ^-ensemble context). In our non-equilibrium 
setup, we work in the path space and condition on the whole trajectory Y = {y(f) : t £ [7j, G]}; 
starting at some time Ti > t\ chosen later. The configuration interval J(t) = [yT-K-i{t),yL+K+i(t)\ 
for the conditional measure is time-dependent, but by rigidity it is quite close to the corresponding 
interval [yitL—K—i){t)> li(L+K+i){t)\ given by the quantiles that remains practically constant owing to 
the removal of the mean velocity. Still, J (f) may wiggle on the rigidity scale N^/N which is much bigger 
than our target scale, 1/iV, the size of the gap, so that we cannot tolerate this imprecision. Furthermore, 
similarly to the basic idea of the local relaxation flow [251 26] we want to achieve universality by showing 
that the measure converges to a (local) reference equilibrium measure. The local Gibbs measures with 
boundary condition y(t) change too quickly to serve as useful reference measures. 

Therefore, in the third step, we define a time-independent local measure, cyri with exterior points 7 *,, 
k £ I c . These exterior points coincide with yk{Ti) for k far away from L while they are given by a 
typical configuration z of an auxiliary quadratic /3-ensemble for k near the boundary of I (with a 
smooth interpolation in between). The auxiliary ensemble is chosen in such a way that the local 
density around A* matches. Using the rigidity bounds for both y and z, we establish that satisfies 
the logarithmic Sobolev inequality (LSI) and the corresponding dynamics approaches to equilibrium 
on a time scale of order K/N. Furthermore, we show that the measure ojt, is rigid by using a general 
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criterion for rigidity of local measures given in Theorem 4.2 of |28| together with the careful choice of the 
auxiliary ensemble. Moreover, we notice that wjy satisfies a level repulsion bound due to Theorem 4.3 
of [2B]. Finally, Theorem 4.1 of [25] implies that the gap statistics of wiy are universal. 

The fourth step is to consider x-i (t) . i £ I, t > Ti, the solution of the local DBM with exterior 
points 7 k, k £ I c , and with initial condition Xi(T\) = ay(Ti). Writing the distribution of x(f) as gt wth 
we derive fast convergence to equilibrium, i.e., for times t>T[ := T\ + K(K /N) the measure gtOJ^ 
is exponentially close to equilibrium in the relative entropy sense. This information can be used to 
transfer rigidity and level repulsion from to g± ujt-i , furthermore it shows that the gap statistics 
of x(f) are the same as those of u>t i: hence are universal. 

The next idea, in the fifth step, is to couple the evolution of x to x by using the same Brownian 
motions in the DBMs. This basic coupling idea first appeared in m in this context; its main advantage 
is that taking the difference of the original DBM and the DBM for x, we see that the difference vector 
v := x — x satisfies a system of ordinary differential equations (ODEs); the stochastic differentials drop 
out. Roughly speaking, these ODEs have the form (see m) 

^ = -(^v)j + Fi , (^v)i := y: Bij(vi - Vj ) + WjVj , (2.28) 

jei 


with time-dependent coefficients B t j , IT) and a “forcing term” F i that all depend on the paths x(f), x(i). 
These coefficients are crudely given by 


_ _ 

lJ N(xi — Xj)(xi — Xj) 


Wl N ^ ( x i ~ Vk){Xi - 7k) ’ 


Vk - 7 k 


p. = _. 

N ^0 ( x i ~ Vk)(Xi ~ 7fc) ' 


The equation (12.281) is very similar to the basic equation studied in f23] but the forcing term is new. 
The key result of [28] is a Holder regularity theory for (12.281) without forcing, under suitable conditions 
on the coefficients. We extend this statement to include the forcing term; here we rely on the finite 
speed of propagation, proved also in |28|. Holder regularity in this context yields that, after some time 
of order K c /N , c ~ 1/100, the discrete derivative Vi+i — iH is much smaller than its naive size 1/N. 
Since 1 — Vi = Xi+± — Xi — (aij+i — Sy), we see that the gaps of x and x coincide to leading order. 
Since the gaps of x were shown to be universal in the previous step, we obtain that the gaps of x(T), 
T :=T[ + K C N~ 1 , are universal. 

There are several technical complications behind this scheme, most importantly we need to regu¬ 
larize the local singularity in the kernel Bij when ay ss Xi±\. In fact, two different regularizations are 
used; the regularization of the dynamics in Section [5.II is borrowed from Section 3.1 of |15j . while the 
regularization of the equilibrium measure cuTy explained at the end of Section [4. 5 1 is similar to the one 
in Section 9.3 of [28] but with a different choice of regularization scale. 


3. Concepts 

In this section we recall essential concepts that will be used in the proof of Theorem 12.11 

3.1. Definition of general /3-ensembles .We first recall the notion of /3-ensembles or log-gases. Let 
N £ N and recall the definition of the set h C in (12.101) . Consider the probability distribution 
on h ' iV ) given by 


4?(dA) = MdA) := J-e“^ w dA, % = U{\) := ^ \v{K) ^ ^ £ log(A, - A. t ) (3.1) 

^ i= 1 Z V l<i<j<N 

where /3 > 0, dA := 1 (A £ F^) dXidX '2 ■ ■ ■ AXn, and Zy = zjfy is a normalization. Here V is a 
IV-independent potential, i.e., a real-valued, sufficiently regular function on R to be specified in each 
case. In the following, we often omit the parameters N and /3 from the notation. We use P w and E w 
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to denote the probability and the expectation with respect to gy. We view gy as a Gibbs measure 
of N particles on R with a logarithmic interaction, where the parameter /3 > 0 may be interpreted as 
the inverse temperature. We refer to the variables (A,) as particles or points and we call the system a 
/3-log-gas or a /3-ensemble. We assume that the potential V is a C 4 function on R such that its second 
derivative is bounded below, i.e., we have 

inf V"{x) > - 2 C V , (3.2) 

x(£R 


for some constant Cy > 0, and we further assume that 

V(x) > (2 + c)log(l + \x\), 


x £ 


(3.3) 


for some c > 0, for large enough |a;|. It is also well known, see, e.g., [12], that under these condi¬ 
tions the measure is normalizable, Zy < oo. Further, the averaged density of the empirical spectral 
measure, g'y' 1 , defined as 




(3.4) 


converges weakly in the limit N —> oo to a continuous function, gy , the equilibrium density, of compact 
support. It is well known that gy satisfies 


V'(x) = ~2 [ 

JR 


Qv(y)dy 

y-x 


x £ supp gy . 


(3.5) 


In fact, equality in (13.51) holds if and only if x £ supp gy. 

Viewing the points A = (Ai) as points or particles on 
cle, 7 *;, under the /3-ensemble gy by 


/ Ik 

gy(x)dx = 

-OO 


k 

N 


we define the quantile of the fc-th parti- 


(3.6) 


For a detailed discussion of general /3-ensemble and the proof of the properties mentioned above we 
refer, e.g., to mna- 

Assume for the moment that the minimizer gy is supported on a single interval [a,b\, and that V 
is “regular” in the sense of [37], i.e., the equilibrium density of V is positive on (a, b) and vanishes like 
a square root at each of the endpoints of [a, b\. From [T3HH] we then have the following rigidity result. 

Proposition 3.1. Let V £ C 4 (R) be a “regular” potential and assume that gy is supported on a single 
interval. Then, for any £ > 0 there are constants Cq,c\ > 0, such that 


' (|A fc - Tfcl > < e 


-c 0 N c i 


1 < k < N, 


(3.7) 


where k := min{fc, N — k + 1}, for N sufficiently large. 


Proposition [3TT] will only be used as an auxiliary result (see Subsection l4.4.2l belowl. since, for most 
potentials of interests in the present paper, the equilibrium density gy is not supported on a single 
interval. The extension of Proposition 13.11 to that settings has not been established. 

Finally, for the Gaussian case, V(x) = ar/2, we write /xg instead of gy, since gc is the equilibrium 
measure for the DBM. More precisely, the Gaussian distribution on h < ' N ' > is given by 


N 


g G ( dA) = —e-^'MA, 
z G 

where Zq = Z^ G is the normalization. 


n G ( A):=El A ?- 1 


2=1 


N 


lo s( A t - A *) > 


(3.8) 


l<i<j<N 


li 





3.2. Dyson Brownian motion. Consider the DBM, A(t), t > 0, on pW of (12.111) with initial 
condition A(0) G f( N ) . Denote by fo p, g, the distribution of A(0 )Q and let f t /ic denote the distribution 
of A(f). Then f t = satisfies the forward equation 


where 


dtft = seu , 


AT 


Jzf = JSf, 


1 


AT 


AT 





(3.9) 


(3.10) 


or in short 2z? = iA — (VHg) • V, with "Hg as in (13.81) . 


3.3. Relative entropy, Bakry-Emery criterion and the logarithmic Sobolev inequality. A 

cornerstone in our proof is the analysis of the relaxation of the dynamics (13.91) . Such an approach was 
first introduced in Section 5.1 of [23]- The presentation here follows |26| . 

Let /i be a probability measure on pW be given by a general Hamiltonian "H: 

d/i(x) = 4e _/3JVH(x) dx, (3.11) 

Zj 

and let L be the generator, symmetric with respect to the measure d/r, defined by the Dirichlet form 

D{f) = DM) = -J fLfdfJi ^ E / (^/) 2d /- i - d i = d *r ( 3 ' 12 ) 

The relative entropy of two absolutely continuous probability measures on R w is given by 


S(is\v) 



dp. 


If dp = f dp, then we use the notation S v (f) := S{f v \ v )- The entropy can be used to control the total 
variation norm via the well-known inequalities 


j\f-i\dp<MsM), p fI '(A)<r'(A) + MsAf), 


(3.13) 


for any immeasurable event A. 

Let ft be the solution to the evolution equation dtft = Lf t , t > 0, with a given initial condition fo- 
Assuming that the Hamiltonian "H satisfies 

V 2 H(A) = HessH(A) >'0, A <E F (A °, (3.14) 

the Bakry-Emery criterion [5] yields the logarithmic Sobolev inequality (LSI) 

£>(/)< 7), f = f 0 eL°°(R N ,d\), (3.15) 

and the exponential relaxation of the entropy and Dirichlet form 

SM) < e- Mt SM 0 ), Dpiy/Tt) < e- 2W D,(t/A), t > 0. 


We assume from now on that H is given by Hg (see (13.81) 1. L is given by 2z? (see (13.101) 1 and that 
the equilibrium measure is the Gaussian one, /i = /iq ■ We then have the convexity inequality 


v,v s h g (x)v) > l iivf + iy; hi-Al > 1 iivi 


Ui - XjY 

1<J v J ' 


v G 


pAT 


(3.16) 


This guarantees that /tg satisfies the LSI with •d = 1/2 and the relaxation time is of order one. 


* Strictly speaking, the distribution of A(0) may not allow a density fo with respect to //q, but for t > 0, A (t) admits 
such a density. Our proofs are not affected by this technicality. 
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3.4. Localized measures. Following [25], we choose K £ [iV ro , iV 1 / 10 ], for some small 0 < vo < 1/10 
and pick L £ \K. N — K\. We denote by / = \L — K, L + K ] a set of /C := 2 K + 1 consecutive indices 
around L. Recall the definition of the set F^ C in (12.101) . For A £ F^ N \ we rename the points as 


A = (Ai, A 2 ,..., A n) = (j/i, • ■ ■, Vl-k-IjXl-Ki • ■ •, xl+k, Vl+k+i , • ■ ■, 2/jv) , (3-17) 


and we call A a configuration (of AT particles or points on the real line). Note that on the right side 
of (13.171) the points retain their original indices and are in increasing order, 

x = (x L _ K ,..., x K+L ) £ F (c) , y = (yi, ■ ■ ■, Vl-k- i, Vl+k+i, ■ ■ ■, Vn) G F {N ~ K) . (3.18) 

We refer to x as the internal points or particles and to y as the external points or particles. In the 
following, we often fix the external points and consider the conditional measures on the internal points: 
Let v be a measure with density on F^ N ^ ■ Then we denote by v y the measure obtained by conditioning 
on y, i.e., for A of the form (13.171) . 


iA(dx) = i/ y (x)dx 


z/(A)dx i/(x, y)dx 
fv( A)dx f v(x, y)dx ’ 


where, with slight abuse of notation, i/(x, y) stands for i/(A). We refer to the fixed external points y as 
boundary conditions of the measure v y . For fixed y £ F (Ar_K 3, all ( Xi ) lie in the configuration interval 


Jy '■= [vl-k-i,Vl+k+i] • (3.19) 

Thus v y is supported on ( J y ) K D F^ , but with a slight abuse of terminology we often say that v y is 
supported on J y . In case v = //i, we define the conditioned density by f y p y = (/p) y . 

For a potential V, we consider the /3-ensemble p.y of m- For K, L and y fixed, we can write p y as 

My( dx ) = ^re^ My(x) dx, R y (x) = ^ ^V y (x^ - log \xj - Xi\ , (3.20) 

v i£l i,j&I 

i<j 

Xi £ J y , with Zy = Zp V a normalization and with the external potential 

V y {x) = V(x) — ^^log \x — yi\ . (3.21) 

igl 


4. Localizing the measures 

4.1. Localization at time T\. Let K £ Nat satisfy (12.271) . with K > |_./V CT J, 0 < vo < 1/10, i.e., 

K W N 5 < N , K>N . 

Recall the constant a > 0 from the assumptions of Theorem 12.11 Let \ G (0,t*7) be a small constant, 
to be chosen later on. Note that N x -C K. Then introduce the intervals of integers 

I := \L — K, L + Kj, I 0 := \L - K 5 - N X I<\ L + K 5 + N x K 4 j , I a := [L — crN, L + aNj , (4.1) 

and we denote by 7 C , Jq, I£ the complements of I, I 0 , I a in Njv- Note that I C / 0 C I c ,. For a 
configuration A £ F^ N \ we introduce x and y as in (13.181) . 

Fix a small £ > 0. Let 

Q\ : = {y G F (JV_/C) : |lfc-7W«)l < ^/^ , Vfc £ J CT } , (4.2) 
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respectively, 


Gt ■■= y e F 


-(N->C) . 


J_ 1 _ f 8s{y)dy 

N y k -x Jye[y L -„ N , VL+aN Y V~ x 


< 


N* 

W 


Vx e [yL-Ks,yL+Ks} /, (4.3) 


where K, = 2K + 1. Note that for each s > 0, we choose the labeling in Q) to be the one of Assump¬ 
tion (2) of Theorem 12. II We then set 

Q s ■■= Q\ n Gs ■ (4.4) 

For any Y := (y(s) £ h ; s £ i, ^ 2 ]} trajectory, we define the conditional measure P Y 

on the X := {x(s) £ jr(JV-JC) . g g [t 1; t 2 ]} trajectories in the usual way. We use P Y to denote 
the conditional measure on the whole X trajectories, while for any fixed time s, we use P>4 S ) for the 
conditional measure (on the x(s) configurations, for any fixed s). We set 


G ■= {Y = {y(s) : s £ [ti,t 2 ]} : 

y (s)£S s , Vs : P Y (|xi(s) — 7^(i)( s )l < N^/N, Vi £ /) > 1 — N~ D , 

and sup max| y k (t) - y fc (s)| < N^y/t - sf , (4.5) 

t,se[ti,t 2 ] fee/c J 


for small £ > 0 and large D > 0. 

Lemma 4.1. Let A (t), t > 0 be the DBM on of (12.111) with fixed initial condition A(0). Let 

Xi(t) = A i(t), i £ I, respectively y k (t ) = A k (t), k fL I, for all t > 0. Under the assumptions of 
Theorem \2. 11 for any (small) f > 0 and any (large) D > 0, we have 

P (y(s) £ Gs, S £ [ii, t 2 ]) > 1 - N~ cD , (4.6) 


and 


P(S) > 1 - N~ cD , (4.7) 

c > 0, for large enough N > No(f,D), where P is with respect the Brownian motions ( Bi ) in (12.111) . 

Proof. Both estimates (14.61) and (14.71) follow directly from the assumptions (l)-(4) in Theorem 12. II □ 

Throughout the rest of this section we will consider the trajectory Y £ Q as fixed. Nevertheless, 
all estimates will be uniform on Q. In particular, all constants only depend on the constants in (14.51) . 
the constants S, e and a of Theorem EH as well as the parameter £ > 0. 

4.2. Regularity of the semicircular flow and removal of mean drift. Consider the DBM, A (t) 
of (12.111) with initial condition A(0) and the semicircular flow p t = Ft\o\- We first study some regularity 
properties of g t for t £ \t\,t 2 \. The following result is proven in Subsection I A. 1.2 1 of the Appendix lAl 

Lemma 4.2. Under the assumptions of Theorem \2.1[ the semicircular flow gt satisfies 

C- 1 < et(E) < C , \d E Qt(E)\ < CN S , (4.8) 

for all E £ [E* — S/2 ,E* + S/2] and all t £ 

Let L £ Njv be as in Theorem EH In particular, we have pti(7L(ti)) > c > 0. Then, we have 

\xi(t) < -jjr, *£/, y {t)&Gt, (4.9) 
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with high probability for N sufficiently large, uniformly in t £ [^i, ^ 2 ], for some labeling £ that will be 
fixed throughout the paper. Recall from (12.91) that the quantiles 7 are determined by 


rik{t) 


Qt(x)dx = — . 


(4.10) 


The evolution of 7 is studied in the Appendix [A] where the following result is proved. 

Lemma 4.3. Under the assumptions of Theorem \2. 1[ the quantiles ( 7 *,) defined through (14.101) . satisfy 

Qt{y)dy 7 


7<(»)(*) = - / 
Jr 




t{i) G la , 


(4.11) 


for all t > 0, where 7 eufit) = gj 7 ^(i)(i). In particular, we have | 7 f(i)| < C. Further, we have 


1 N 

T’iWW = 


7 


N f-' 7fc(t) - Uc{i){t) 2 


k—1 


°( 


N 6 \ 

“iv J ’ 


£(*) G la , 


uniformly in t £ \T\ , f 2 ]. Moreover, we have the estimates 

It/«>(*) - 7/(L)(i)l < CA- 1+5 |* - L|, |7/(i)W - < CN*(t - Ti), 


(4.12) 


(4.13) 


uniformly in t £ [Ti,£ 2 ] and £(i) £ I a . 

Equation (14.111) shows that the points (7 e(i)(t)) approximately satisfy a gradient flow evolution of 
particles with quadratic confinement and interacting via the mean field potential log \ x — y\. 

Lemma S3] is proved in Subsection I A. 1 .2 1 of Appendix [A] Let us briefly mention how the constant E 
in Assumption (1) and the constant a in Assumption (2) can be related. For given E > 0, we can 
choose a > 0 such that 7 ^) (t± ), for any i £ I a = \L — aN, L + ulVj, all lie inside [A 1 * — E/4, A* + E/4]. 
Then we know from Lemma FCil that 7 t(i)(t) £ [A* — E/2, A* + E/2] for all t £ [fi,f 2 ]. By Lemma 14721 
we have control over g t on [A* — E/2, A* + E/2]. 

For simplicity of notation, we henceforth drop the labeling £ and simply write, with some abuse of 
notation, 7 fit) = 7 From (14.131) and (14.91) . we conclude that 


\X L (t)~ Ai(Ti)| <v L {t-T{) + o{t 


T\) + O 



t £ p~i,£ 2 ], 


with high probability, where we have set v /, := 71 , (Ti). We denote by A (t), t>T\ the process obtained 
from A {t) by setting 


A(£) := A(i) - Ul(£ - Tf), t > Ti. 


(4.14) 


Thus A(f) satisfies the SDE 


d\(t) 



v L d t 


1 1 ^ Aj(£) + wl(£ — Ti) 

N ~?i *j(t) - *»(*) 2 


i£N N , /? > 1, 


for t > T\. In the following we write Xi = A*, i £ I, respectively y k = A&, k /, so that Y = 
|y(t) £ : t £ [Ti,t 2 ]}. Having shifted the original process (A(£)) as in (14.141) . we also shift the 

distribution g t and the quantiles 7 , for t > Ti, accordingly: 


g t (x) := q t {x + v L (t-Ti)), 7iW : = 7»W - v L (t - T{) , tG[Ti,i 2 ], 


x £ R, i £ Nat. In a similar way, we introduce the events Q s , Q s , Q s and Q by replacing the quantities 
without bars with bars in (14.21) , (14.31) , (14.41) and (14.51) . 
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4.3. The reference points jj. Once Y £ Q, thus also Y £ Q. is fixed, we introduce time-independent 
“reference points”, 7 = ( 7 *,), as follows: For k £ N^v, let 

L - K 5 - N X K 4 <k<L-K 5 , 

L + K ' 5 < k < L + K 5 + N X K 4 , 

L — K 5 < k < L + K 5 , 

) 

with x > 0 as in (EH, i.e., ik is a linearly mollified cutoff of the indicator function 1 (k £ I a ). Set 

Ik ■= Ik Zk + (1 - Lk) yk{Ti) , k £ I c , (4-15) 

where the external points z = (zk) £ in (14.151) will be chosen in Subsection 14.4.21 below. Note 

that 7 k = Zk, for|L — k\ < K 5 ; 7 t = yk{T\) = y k [T\), for \L — k\ > I\ 5 + N X K 4 . Thus the sequence 7 

smoothly interpolates between the external points y(Ti) from the DBM and z. The external points z 
are constructed from an appropriate /3-ensemble whose equilibrium density has a single interval support. 
This will guarantee rigidity; in particular, 


L’k * — 


■ k-(L-K 5 -N x K 4 ) 

NxK 4 11 

L+K 5 +N x K 4 -k ^ 


NXK 4 


if 

else 


Ufc+i ~fc| < C ^ 2/3 jzi/z ’ 

with /c = minjfc, TV — fc + 1}, for all k £ [1, IV — 1]; see (14.301) below. 

Anticipating the precise choice of z, we mention that they are chosen such that 


(4.16) 


ZL—K—i = , zl+k+i = 2/L+if+i(7i). (4-17) 

In fact, this choice will assure that the configuration interval of the localized measures, both with y(Ti) 
and with 7 as external points, will have the same (and time-independent) support 


Jy(T{) — Ji. — [z~, z+], 

where z_ = zl-k- 1 , 2 + = Zl+k+i- We next estimate the size of the interval J z . 
Lemma 4.4. Let J z be as in (|4.18l) and assume that K satisfies (12.271) . Then, we have 

IC 


(4.18) 


\J zl = 


Ng Tl {z 0 ) 




(4.19) 


on QT it where z Q := {z l +k+ 1 - z L -k- i)/2 . 

Proof. We mainly follow the proof of Lemma 4.5 in [5S]. First, we write, by (14.181) . 

\Jz\=z + — Z- = yL+K+i(Ti) - y L _ K -i(Ti ) = 7 l+k+i{T \) — + 0(TV~ 1+ ^), 

where we used that y(Ti ) G Gk- Next, we note that by Assumption (1) of Theorem 12.11 we have 
BtA x ) = £Ti (z) + O ( N 5 \x — z\). Thus from (12.91) . 

/•7l + k + i(Ti) 

IC = N g Tl (x) dx = Ng Tl (y)| J z | + 0(TV 1+<5 | J z | 2 ) -)- O(TV^), 

J JL-K-i(Ti ) 

where we used that qt x ~ 1. Since I\ N s <C TV, by (12.271) . we get (14.191) . □ 
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4.4. Localizing the DBM and the reference measure Having Y £ Q fixed, we consider 

the DBM on the x-variables given by the stochastic differential equation (SDE) 


d Xiit) 



v L d t 


- Y- _ 1 _ 

N Xj(t)-Xi{t) 

j¥=i 


d t 


1 

N 


E 


i 


d t — 


Xj(t) + vl (t — Ti) 


(4.20) 


i £ I, t > Ti, with a collection of independent standard Brownian motions. We let P Y denote 

the associated path space measure. 

For t > Ti, we define an approximate coupled dynamics, x(t) by letting 


d Xi{t) 



v L d t 


i i dt L L d t 

N -Xi(t) 1 N “ 7 fc — Xi(t) 


E^dt, (4.21) 


i € I, with initial condition x(Ti) = x(Ti). The corresponding path space measure is denoted 
by P Y . Going from (14.201) to (14.211) we replaced the time-dependent external points y (t) by the 
time-independent reference points 7 and we neglected the drift term 17 (t — Ti)di/2. Note that the 
Brownian motions in (14.201) and (14.211) are the same. 

We define the local “reference” measure 


wii(x) dx := J-e-f* NHT iW dx, ft Tl (x) := £ V^x*) - ^ ]T lo g(^ " **) 

Tl iei i,jei 

i<j 


(4.22) 


where the external potential is given by 


V~<(x) := \v{x) - ^^]log|x- 7 fe | 


k&i 


V(x) = — + 2v L x. 


(4.23) 


The subscript T\ in indicates that the external points 7 in the construction of this measure were 
obtained in (14.151) by matching the external points y(Ti) of the original DBM at time T\. Note that 
this measure as well as the measure P Y are supported on the fixed configuration interval 


*Jy(Ti) = [yL-K-l{Tl), yL+K+l{Tl)\ . 

The measure wtj is the equilibrium measure of the SDE (14.211) . 

We write the distribution of x(t) as gtOJr i (for t > Ti). Since they are supported on the same 
configuration interval, the measures gt Wti (for t > T\) and wjy are both absolutely continuous with 
respect to the Lebesgue measure, hence also to each other. 

4.4.1. Entropy bound .In this section, we compare the measures wji and g t u;Ti for t>T\. We show 
that the process (x(f)) equilibrates on a time scale ~ K/N , i.e., the local statistics of gtuiT! an d wji 
are very close beyond times t > T[ := Ti + K(K/N ), with t < £ 2 - 

Since is supported on an interval of size 0{K/N ) (see Lemma T4.4L the Hessian of its Hamil¬ 
tonian "Hti from (14.221) satisfies 


njH 

n Tl 


, , . 1 V- 

(x) > min — 


Lilli! - 7 — - 

N {X 


1 . 1 

-——-77 > mm — 

- 7 k) N 




> 


k : K<\k—L\<K 5 


(Xi - Zk) 2 I< 


77 - ( 4 - 24 ) 


for all x e H F (,c \ where we used (14.161) . Thus, recalling the discussion in Section 1X51 

satisfies the logarithmic Sobolev inequality 


CK 

S UTl (f)<—D, u 


(V7); 


(4.25) 
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c.f., (13.151) . Further, for t > T\ + i K(K/N ) the process (x(t)) has become absolutely continuous with 
respect to Lebesgue measure, and one can easily prove that 

S(g t wr, \oj Ti )<N C , t>T x + \k{K/N) ; (4.26) 

for some large C; see, e.g., Lemma 4.7 in |21| . Therefore, running the Bakry-Emery argument of 
Subsection 13.31 from time T\ + 1 K(K/N ) to time T( = T\ + K(K/N) and using the initial entropy 
estimate (14.261) . we immediately get the following result. 

Lemma 4.5. For any t > T[ = T\ + K(K/N), we have 

Du, Tl {V9l) + S UTi (g t ) < e~ cK , t G [T[, t 2 ] , (4.27) 

for some c > 0. In particular, the statistics o/x(f) for any t G [T[,t 2 ] are the same as the statistics of 
the local equilibrium measure as follows from 


/ 


0(g t 


1 ) d(jUTi 


<\\0\\ xy /2sZJgi)<Ce- cK /\ 


(4.28) 


for any bounded observable O. 


4.4.2. Construction of an auxiliary /3-ensemble.We now turn to the choice of the reference points z 
introduced first at the beginning of Subsection 14.31 We construct a global /3-ensemble, /r AUX , with 
potential V1 UX and equilibrium density p AUX such that it has a single interval support and such that the 
density matches with qt, at ')l{T\). The main properties of /r AUX are summarized in the next lemma. 

Lemma 4.6. There exists a /3-ensemble /r AUX = /r lux, with quadratic potential b AUX and equilibrium 
density g Am , and a set of external configurations z G /r(W— k:)^ w ^ Zl k i = y L _ K _ 1 (T\), zl+k+i = 
yL + K+i(Ti), such that the following holds for N sufficiently large. 

(1) The limiting equilibrium density g Am of g Am is a shifted semicircle law with finite variance 
satisfying, for any £ > 0, 

£aux0/) = Qr^y) + 0(N*/K) + 0(N s \y - z L -k-i\) , 2 / G M. (4.29) 

(2) The external points z satisfy, for any £ > 0, 

l* fc - 7 Aux, fe |<C'^p kel c , (4.30) 

where (7 AUX ,fc) are the quantiles of the equilibrium density g AUX , i.e., Jf AUX,k dg AVX = k/N, and 
k = min {k,N — k + 1}. In particular, since 4^ ux is “regular”, the rigidity estimate (14.1611 holds. 

(3) The localized measure /r AUX satisfies, for any £ > 0, 

P ^™(|Xi - oti\ > N*/N,Vi G I) < CW 5 — , (4.31) 

where ( at ) are K. equidistant points in J z = [zl-k-i, Zl+k+i] = J y (Ti)> he., 


Oti = z Q 


i — L 
2K+1 


Zo = —{zl-K- 1 + ZL+K+ 1 ) ■ 


(4.32) 


Proof. The proof is split into three steps. Step 1: We introduce the quadratic potential Vq/^ (x) := 
x 2 /2<i 2 , with some q > 0, and consider the /9-ensemble, gG(s)i with Hamiltonian 


n G{s) ■= 


1 

2 


1 V— v 

— J2 log |Aj — Ai|. 

i,j =1 

i<j 
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It is easy to check that the limiting equilibrium density, Qg(s)> of /^g(c) satisfies Qgu)(x) = Qsc(x/<;)/<;, 
with Q ac {x) = f \/(4 — x 2 )+ the standard semicircle law. Similarly, the quantiles, (7 g( 0,*)> °f 0G(?) 
satisfy 7 g(^),i = ?7 sc,i, where 7 SCj i denote the quantiles with respect the standard semicircle law, 
i.e., dg sc = i/N. Thus £>g(.;)(7g«m) = Qsc{lsc,i) /?• In particular, we can fix ? such that 

0t 1 (7l( t i)) = 0G(?)(7G(?)j)> *-e., we set 


Qsc(Sfsc,L ) 

£'Ti(7i(T , i)) 

We next choose boundary conditions y with the following properties: (1) For any £ > 0, 

li/fc - 7sc,fc| > N*/N, Vfc G I c , (4.33) 

for AT > ATo(£), (ie., y are rigid in the sense of sense of (13.711 with V = Vq(?)); (2) for any £ > 0, there 
are c ' 0 , c' x >0 such that 


P^wfl Xi - a*| > N*/N ,\/i e I) < e~ c '° N '’' 1 , (4.34) 

where cL are the IC equidistant points in the configuration interval J y = [i/L-ar-i, Vl+k+i]- The 
precise choice of y is unimportant for our argument, as long as y satisfy (14.331) and H4.34[l . That we 
can choose a y such that (14.331) and (14.341) are satisfied follows from Proposition 13.II and an application 
of Markov’s inequality. 

Step 2: The length of the configuration intervals J y (r 1 ) and J y may differ slightly. Using the 
scale invariance of the Gaussian measure, we now adjust ^ and y to guarantee that the lengths of the 
configuration intervals agree: Following the proof of Lemma U~T1 or the proof of Lemma 4.5 in |28j . we 
get from the rigidity estimates for Pg(c) that 


\Jy\ = \vl-K-i - VL+K + 1| = 


K. 


N QG(s) (7g (?)j) 


and from Lemma 14.41 that 

l*^y(7i)I = \ul-k-i{Ti) - Vl+k+i(Ti)\ = 


IC 


0(N~ 1 N ^), 


+ 0(N~ 1 N i ). 


NqtAil) 

Using that Qt x (7l) = f?G(?)(7Gby our choice of we hence conclude that 


s := 


= 1 + 0{N^K~ 1 ). 


'y(Ti)l 


(4.35) 


Setting z := y/s we have | J^\ = |J y ( Tl )| and 

P^g) (| Xi- &i/s\ >N*/N,VieI) = P m °g'> (\xi-&i\>8Nt/N,VieI), (4.36) 

where we have set <;' : = Using the rigidity of y we get, similarly to (14.241) . that V^ y (x) > 
for all x € (Jy) K DF K - Thus the logarithmic Sobolev inequality 

CK 


^GU'jl^GU)) - jy 


with the local Dirichlet form 


[= mZj ( a (H^) w) 
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holds. A straightforward calculation together with (14.3511 then shows that 


iei J 


2 - N 2 ^ K 2 

d n* < C- 


‘GW - 


N 2 


Thus, using first (14.361) and then the entropy inequality (13.131) . we get 

r4 w (\ Xi -ai/s\ > N*/N,Vi £ I) = P"«fr'> fla* -a ( | > sN^/N,Vi £ I ) 


< (|*i - &i\ > sNt/N,\/i £ /) + v /25(^ (? , ) |^ w ) 

(4.37) 


<Ce-"‘ + C NtK 


N ’ 


where we used (14.341) (with an additional factor s) to get the last line. 

Step 3: Finally, we achieve that J z = «/ y ( Tl ) by a simple shift in the energy: we replace RgGoC 2 -) 
by Vq(c'){x — b), b := yi,-K-i(Ti) — pl-k-i, x £ K. We now choose /z AUX as the Gaussian measure 
defined by the potential Vg(«')( - — b) and we set Zi := Zi — b , for * S [1, 77]. With these choices, (14.37[) 
asserts that 


P^ ux (|Xj 


-Q.i\ >N*/N,Vi€l) <C 


N^K 

N 


(4.38) 


where are the K. equidistant points in the interval J z = J y (Ti) = [yL-K-i(Ti),yL+K+i(Ti)\. 

In sum, we have established the following. We consider the ^-ensemble p Am with quadratic po¬ 
tential, whose equilibrium density g AUX is a semicircle law with radius \/ 2 <j / which is centered at b. 
Taylor expanding the densities qt x and p AUX around PL-K-iiTf) and recalling (14.351) as well as (14.81) . 
we obtain (14.291) . This proves statement (1) of Lemma [4.61 The points z = (zf) are rigid as follows 
from (14.3311 and the choices z, = yi/s, Zi = i) — b, i £ Nn- This immediately implies statement (2) 
of Lemma 14.61 Finally, the rigidity statement (3) of Lemma 14.61 for the localized measure /i AUX was 
obtained in (14.381) . This concludes the proof of Lemma 14.61 □ 

We conclude this subsection with a straightforward technical result that will be used in the next 
section. Recall the definition of the interval of integers /, Iq and I„ in (14.11) . 

Corollary 4.7. Let z be as in Lemma \4-b] and let 7 be defined as in (14.151) . Then, for any £ > 0, 

\lk - 7 *m)| < C^r + C ■ 1 (k £ Jo) (^ N ^ k K 7l| + N s \% - 7 l| 2 ) , k€l„\I, (4.39) 

for N sufficiently large. Moreover, we have, for any £ > 0, 

7 k - 7k -1 < N*/N , k £ I a \I , (4.40) 

for N sufficiently large. 

Proof. Recall that 7 *, = ykiTf) for k ^ I 0 . Since y(Ti we immediately get 


|7 k ~ Jk(Ti)\ < C 


Nt 

~N~ 


(4.41) 


for k £ . Next, assume first that A' + 1 < L — A: < K 5 . Then we have 7 *, = Zk, and we can write 

O 


fyK- L - i(Tx) \L-K-l-k\ 

0Aux(?/)dy = — 

Ilk 


N 


7V«\ 

W 


20 






























where we used Uk-l-i{Ti) = Zk-l- i, the rigidity estimate in (14.3011 and the fact that (7 AUX ,fc) are 
the quantiles of p AUX . Using (14.291) . we hence can write 


rVK-L-i(Ti) 


Qn (y) d y 


\L — K — 1 — k\ 0 (N*\ 

N \N J 

+ o ^ N ^-yL-K-l(Tl)\ ^ + Q _ y L _ K _^ Tl )\ 2 ) . 


On the other hand, since y(Tj) G Q^ ± , i.e., \yL-K-x{Tf) — 'Yl-k-i(Ti)\ < CN^/N, and using 
that 7 fc(Ti) are the quantiles with respect to qt x , we have 


rVK-L-i(Ti) 


hk(T i) 


QT X (y) d y 


\L — K — 1 — k\ 
N 



Comparing these last two equations and using the lower bound on the density pji , we conclude that 


\z k - 7 fcCT 1 )| <C^ + C N \ 7l| + CN S | 7 fc - 7 l | 2 , (4.42) 

for k such that K + 1 < L — k < K 5 . Here, we also used that 7 l — t/£_^_i(Ti) < CK/N. The same 
argument applies to the case K + 1 < k — L < K 5 . 

It remains to consider the transition regime A' 5 < \L — k\ < K 5 + N X K 4 . Using the definition of 7 
in (14.1511 . we can estimate 


|7 k - 7fe( r i)l < t'k\zk - 7fc(^i)l + (1 - tfe)|j/fc(Ti) - 7 fc(Ti)| 


< C 


N* „ 1 7 ^ — 7x,| 


N 


C 


K 


+ CN 5 |7 fe -7L| S 


for such fc, where we used (I4.41|) . (|4.42l) and the rigidity of y(Ti) G • This proves (14.391) . 
The estimate (14.401) follows directly from the rigidity of y(Ti), z and (14.391) . 


□ 


4.5. Three measures and their properties. Having Y G Q fixed and having constructed the 
external points z, we have, up to this point, introduced three distinct measures on the internal particles: 

(1) tOTi is given by an explicit formula in (14.221) . It is a local /3-ensemble on J z which we refer to as 
the “reference” measure. 

(2) gtuiTi is the distribution of x(£) from the dynamics (14.211) on J z . 

(3) P y W is the measure of the x(t) dynamics (14.201) at time t. it is also the conditional measure P Y of 
the original measure P, conditioned on the Y-trajectory at time t>t\. This measure is also on K. 
particles, but now the configuration interval is time-dependent Jy(t) := L-K-i(t)>VL+K+i(t)]- 

In the remainder of this subsection, we establish rigidity for the measures u>t i and 7 * ujt x • 

Definition 4.8. We say that the measure v (on IC-point configurations labeled with I, |/| = 1C, in a 
fixed interval J) is rigid with exponent f if 

i/flxi - a 4 | > N*/N, Vi G I) < Ce~ cN( , (4.43) 

where ai are the K. equidistant points in J and where c > 0. The path-space measure Q for times 
[Ti, £ 2 ] on the same configuration interval J is rigid with exponent f if 

Q( sup |xi(a) - «i| > N*/N, \/i G /) < Ce~ cNi . (4.44) 

sG[T! ,i 2 ] 

Note that if for all t the fixed-time marginals Q t of a space time measure Q satisfy rigidity, then Q 
satisfies rigidity (since the trajectories typically have some mild continuity; see Section 9.3 of |28|). 
We will establish the following main technical input. Recall that T[ = Tf + K{K/N). 
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Proposition 4.9. Let £ > 0 be sufficiently small and let I\ satisfy (12.271) . Then, for any Y £ Q and 
any t £ pi,£ 2 ] the following holds. 

( 1 ) wti fi.e., the local “reference’’ measure) is rigid with exponent £ and satisfies 


maxE“ T i 

i£l 


1 

[N\xi - Xi±i\) p 


< C v 


(4.45) 


(with x i± i(t) = yL±{K+i)(Ti) if i = L ± K), for any p < 2. 

(2) g t ujt x ( i-e the time marginals of the process X = {x(s) : s £ \T[, tf\}) is rigid with exponent 
moreover, the whole process {x(s) : s £ [X^,^]} with measure P Y is rigid with exponents £. 
Furthermore, 


maxE 9 * WTl 
id 


1 

[N\xi(t) - x l± 1 (t.)\]P 


< C p , 


(4.46) 


(with Xi±i(t) = Vl±(k+i)(Ti) if i = L ± L), for any p < 2 and t>T[. 


To simplify the exposition, we split the proof of Proposition 14.91 according to its statements. 

4.5.1. Proof of statement (1) of Proposition^) 7/| We start with the rigidity of the reference measure wji • 
For notational simplicity, we write in the following 


Ik = 7<?(fc)( T i) > k€l a , 

where the labeling i is chosen according to (14.21) . 

Proof of rigidity of ojt, • We first recall the following general result of Theorem 4.2 (see also the remark 
after Lemma 4.5) of }28| . For any local equilibrium measure p y on 1C points with potential V y on an 
interval J of size |«7| ~ K/N rigidity (with exponent £ > 0) in the sense of Definition 14.81 holds if the 
following two conditions are satisfied: 


and 


(^ y )'(z) 


0 ( 3 /o) log 


d +(x) 

d- (x) 


+ 


7V C « 

Nd(x) 


(4.47) 


E^ 


< N^/N, 


Vie I, 


(4.48) 


where y a is the midpoint of the interval J, d(x) is the distance of x to the boundary of J , 
d-(x) := d(x) + g(y a )N i /N , d + (x) := max{p - y L -K-i\, \x - Vl+k+ i|} + 0 (y o )N i /N , (4.49) 


and at are the K. equidistant points in J. 

We now apply this result with the choices y = 7 and J = J z = Jy(Ti) to the reference measure ujtx- 
The condition (14.471) will follow from the global condition (14.521) below and from the fact that the 
reference points 7 are rigid in the sense of Corollary 14.71 The details are as follows. 

To check condition (14.471) . we introduce the supplemental potential V 1 by setting 


V^(x) := i x 2 + 2v l x - 


Y log p- 7*1 


We then have 

(VY(X) (vY(x) 


k:\k-L\>K+Ni 


1 1 N* 

~ N 1 7 * — x\ ~ Nd(x ) ’ 

k-.K<\k-L\<K+Ni 1 ' 1 W 


oc G J z • 


(4.50) 


(4.51) 


To control V 1 , we can follow Appendix A of [25]. First, recall from Lemma POl that vl satisfies 
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7 l = 1l{T\) • 


(4.52) 


VL=lL = 


QtAv) dy _ 7 l 
y~lL 2 


Thus, 


, X f g Tl {y)dy 

VL + „ + / 

_ 

.7£ . f QT,(y)dy 
VL+ „ + / 

2 Jr y-x 


1 

M 


+ 0{\Rem Tl {'yL) - Remr^ar)!) + 0(\x - 7l|) 


AT 5 AT K 

< C——— + C — , 
N N’ 


x e J z 


(4.53) 


where we used (14.521) . To bound the second and third term on the right we used Assumptions ( 1 ) of 
Theorem 12.11 and the estimate on |J Z | in (14.191) . We may thus split 

; K 


(v^y = fii + n 2 + n 3 + o ( n 6 - 


with 


fli(a;) := - f 

^(z) := - f 

J 7 

n 3 (a;) := - f 

J -y 


iL+K + Nt. (y)d y 

y * y - x 

'l-k-n£ * 


^ l-k-nZ 


71 

IN 


Qt, ( y)dy , 1 
y-x 

Qt, (y)dy 


L-K-N* 


E 


iv ~ 7 fc - * 


1 N 

+ v E 


1 


A^ , 7fc — x 

fc=i+if+jve 


'Tt+K+ws y x 

To estimate ffi we use that qt,{ y) = Qt,(x) + 0 (A^ 5 | 2 / — x|) (c./., (|4.8p and (14.191) 1 to get 

x € J z . (4.54) 


ni(z) = -er, ( 2 / 0 ) log + 0{N S K/N) + O 


d_ (x) 


/ jv 2 « 

\N' 2 d{x) 


To obtain the third line, we used 

7 L+K+N* — X = (7 l+K+NS ~ iL+K+l) + {lL+K + 1 ~ x) = d + {x) + 0(N^N X ) , 

respectively x — 7 t+k+n* = d-(x) + 0(N^N~ 1 ), where we used the definition of 7 in (14.15P and the 
definition of d± in (14.491) . 

We next estimate fi 2 (^3 is estimated in the very same way): We split 


L-K-N* 


M—l 


- y ~ 1 =-e 

N 87 - <r /V 


1 


N ti 'Tfc "* ^ IS » fc ( Tl )" * ^ ^ T'fc-* 


. L-K-N* 1 

E 


with M = L — aN, such that we can estimate on one hand 

1 r 7 "- 1 PTi (j/)dz/ 


1 M—l 

- E 

/V E^ 


N y k ^)~ x 


h 1 y - 1 

since y(Xi) £ Qt i) c./., (14.31) . On the other hand, we estimate 
L-K-N* „ L-K-N* 


/AT? 
° ( A^ 


(4.55) 


(4.56) 


1 

IV 


E E 


0T, (y)dy 


k—M 


Ik- X 


k=M J ^k-i 


Ik - X 


fTL-K-Nt 


Qt, {y)dy 
y-x 


O 


1 L-K-N* 

E 


lk \ik — y\or 1 (y)dy 

k=M J Tk-i 


(■ y - x? 


(4.57) 
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Using Corollary 14.71 and recalling the definition of I from (14.11) . we can bound the remainder in the 
above equation as 


L-K-N* 

E 


k=M •'T'fc-i 


7fc \lk-y\0rMdy f 1 l-k-ns d y 

(y-x ) 2 ~ N 


C 


7M-1 (y-x ) 2 

Pl-k-nS (Nt/K)(y-x) + N s (y-x ) 2 


'^L-K^-NXK 4-1 


(■y - xf 


d y 


N* 1V« 

iC mr) +c i<^ 


X - 'JL-KS-NXK 4 -l 


K 5 

+ CN*—. (4.58) 


X — 7 L-K-Ni 

Thus, using that d(x) < |«/ y (Ti)| < CK/N, he., K > cNd(x ), and that Ji satisfies (12.271) . we have 

L-K-N* 


E 

k—M 


lk 17fc - y|£’T 1 (j/)dy < ^ iV 25 


(y - xf 


Nd(x) ’ 


where we bounded the logarithmic term on the right side of (|4.58l) by N*. Hence, combining this last 
estimate with (I4.56|) . we find 


|fi 2 (z)| < C 


7V 2 « 

Nd(x) 


Nt N 2 * 

N s ~ Nd(x) 


(4.59) 


where we used that K > cNd(x) and that K satisfies (12.271) . The same bounds holds for U 3 . 
Combining (14.591) . (14.541) and (14.511) . we get (14.471) for (V^)’ (with c = 2). 

It remains to check (14.481) with the external points y = 7 , i.e., (E^^Xi — oti\ < N c ^/N, i £ I. 
First, we note that from (14.381) we have — on | < CN^/N. Then, using the logarithmic Sobolev 

inequality (14.251) . we bound the relative entropy 


I^tJ < c^Y^ |a i e-^S 1 [ v-(x 4 )-^(xO]| 


N N 


i&I 


< CKW™ Y 


iei 


\vL x(*i) -\xi-v L - ^ Y 


k : \k—L\>K 5 


1 


1 


Sik-Xi z k -Xi_ 


(4.60) 


Note that when k is close to the interval I in the summation above, i.e., when \k — L\ < I\ 5 , then the 
corresponding terms exactly cancel by the choice of 7 in (14.151) . 

To bound the right side of (14.601) . we first recall that we have from (13.51) the equilibrium relation 

Kux(x) = -2 [ gAUX ^ Ay, x £ supp p AUX , (4.61) 

Jr y-x 

for the auxiliary /3-ensemble p AUX . We denote by ( 7 AU x,i)£Li the quantiles of the measure p AUX and let 
7 aux,o = Oaux, 7aux,jv = b Am , where a AU X , 6 AU X are the endpoints of the support of p AUX . 

We then bound the summation in (14.601) for all k < L — K 5 as follows (the case k > L + K 5 is 
treated in the very same way), 



Q Aux(y) Ay 

Zk - Xi 


QAvx(y)Ay \ 

y-xi ) 


L-K” 


<CY 

t, '7AUX,/e — 1 

^5 


L—K 


E 


7A U x,fc \y_ z k \g Am (y)dy 

-1 

N* 


(y - x^ 2 


fe =1 


AT2/3fcl/3 


7AUX,fc &Aux(y)Ay <c Ni 


Taux 


t k~x (y - Xi) 2 k 5 


(4.62) 
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for all i £ I, where we used the rigidity of z (see (14.301) 1 and that p AUX vanishes like a square root at 
the endpoints a AUX , b A ux of its support (recall that p AUX is a rescaled and re-centered semicircle). 

On the other hand, reasoning exactly as in (14.541) . we find that 


pAVX,L + K 5 
^aux,L-K 5 


gAux(y) _ P™x- L + K5 g AVX (xj) + Q(\y 
y~Xi Jry y~Xi 

I aux,L — K° 


■E d y = 0(Jv- 4 ) + 0(K 5 /N ), 


for i £ I. We therefore get, combining (14.611) . (14.621) and (14.631) . 


(4.63) 


\ v LAxi) + ^ 


N z —' Zb — Xj 

k : \L — k\>K 5 * 


< 


c_ 

K 4 ’ 


i G I. 


Second, using the definition 7 in (14.151) . we obtain similarly to (14.551) and (|4.56D . 


(4.64) 


1 

N 


L-K 5 

E 


k= 1 



pL-K 5 
'll 


(y) dy 

y-Xi 


< 


1 

iv 


L — K 5 

E 


k—M 



pL-K 5 
'7m 


g?i fa) dy 

y-Xi 


' C 1W : 


(4.65) 


with M = L — crN. The first term on the right side of (14.651) can be controlled, similarly to (14.571) 
and (14.581) . as 


1 

N 


L-K 6 


E 

k—M 





’ 1M 


QT! (y) d y 


y-Xi 


Nt r /V 5 

< c —+CN 5 — 
~ K 5 N 


C 


AT«+x 


where we used — x,| ~ K 5 /N and the assumption on K in (12.271) . A similar estimate holds for 

the summations over \L + A' 5 ,7VJ. Further, repeating the arguments of (14.541) . we get 


h l+k 5 qtAv) 

IlL-KS y~ Xi 


d y = 0(K~ 4 ) + O 



Thus, combining the last two estimates and recalling (14.531) as well as (12.271) we find 


1 

-Xi + v L + 


1 

N 


E 

\k—L\>K 5 



<° + c 

- K 4 


Nt+x 


(4.66) 


Plugging (1T661) and (14641) into (l4A0l) we get S(jj z auk \oj T i ) < CN 2 ^ +2x K~ 2 , which finally leads, in 
combination with (14.381) , to 


P“ Tl (| Xi - ai \ > N^N- 1 , V* G l) < P"™(| n - «i| > N^N~ l ,Vt G /) + ^/2S(p Am \u Tl ) 




N I\ 

where we used (13.131) . Together with the a priori bound \xi — Ui\ < C(K/N), this implies 


(4.67) 


N* K N i+X N i+X 

|E^ Xi - ai \<C— + C- — < C— , i el. (4.68) 

Thus, choosing, e.g., x = £, we get the bound (14.481) for the measure 

Applying Theorem 4.2 of [28] as mentioned at the beginning of the proof, we see that the mea¬ 
sure ojt-l satisfies rigidity with exponent £. □ 

The level repulsion estimate (|4.45l) in statement (2) of Proposition 14.91 is proved using the explicit 
Vandermonde structure of wji • The proof is essentially identical to the proof of Theorem 4.3 in [28] 
given Section 7.2 of |28]. We therefore leave the details aside. 
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4.5.2. Proof of statement (2) of Proposition \ 4■ 9\ The rigidity for g t u>t x , with fixed t > T[ in the sense 
of Definition 14.81 immediately follows from the rigidity for lot x and the entropy estimate (14.271) . Using 
the stochastic continuity of (x(f)) and the rigidity of gt^>T x a sufficiently large set of discrete times, 
we can conclude that P Y itself is rigid; see Section 9.3 of [2S] for details. 

It remains to prove the level repulsion for g t uj given in (14.461) . 

Proof of (14.461) . The level repulsion bound (14.461) follows from (14.451) and the entropy bound (14.271) . 
More precisely, we have to introduce w^*, an £*-regularization in the u>t x measure in the same way 
as in Section 9.3 of [28j. The parameter e* = e~ K will be chosen tiny with a small c > 0. This 
regularization modifies the interaction terms in (14.2111 and in the Hamiltonian T-Lt x ■ In the latter the 
log becomes log £< defined 

log £i (x) := l(x > £*) log(x) + l(x < £*){ log£* + X ^ £ * - ^-(x - £*) 2 } • (4.69) 

This has the property that d 2 log £ (x) is the same, —x -2 , as before if x > £*, but it remains bounded 
by e~ 2 for all x. The Hamiltonian is still convex. The support of the measure is not J z but the 
whole K., but it is still overwhelmingly supported on J z . In particular, ujt x and are close in entropy 

sense, see (9.57) from [28] . 

S(w Tl |<4;) <CK c e\. (4.70) 

As a consequence, by the entropy inequality (13.1311 we may transfer the rigidity bounds from the 
measure wji to the measure i.e., we have 

P^i (| Xi - ai| > N*/N) < e~ K °. (4.71) 


Similar modifications occur in the SDE (14.211) : the (x x — Xj) -1 and also the (xi — 7fc) _1 terms get 
regularized to 

(Xi - Xj)~ < 1 := d x log £< _ ( Xi -xj), 

and they will be uniformly bounded by £“• Now we can prove (14.461) with the regularization, since 
we can use the entropy inequality (13.131) to get 


E 9t 


[N\xi - z l+ i| e J p 


< E“ T i 


[N\xi - x i+ 1 | 


— +e~^2S UJ ' Tl (g t )<C p KZ, (4.72) 


Here in estimating the first term we used that the level repulsion bounds hold for the regularized 
measure 

P aJ ' r i (xi +1 — Xi < s/N) < CK^s 2 , s > iv^£*, 

see (9.58) of [2B], i.e., we have 


7 —---— < C p Rt. 

[N\xi - x i+ i| £ J p 


(4.73) 


The exponential smallness of the entropy S' ( g t ) is proven exactly the same way as the proof of (14.271) . 

Tl 

since the Bakry-Emery type convexity argument remains valid for the equilibrium measure as well. 
This exponential smallness wins over ef p if the constant c in the definition of £» = exp(— K c ) is small. 

Since the only purpose of this regularization is to prove (14.461) , we will not carry the £* superscript 
throughout the proof, i.e., we continue to write ujtx everywhere, although we really mean As we 
have seen, the key input information on u>t x for our whole analysis, the rigidity (14.711) . holds for the 
regularized measure. The other input, the level repulsion in the form (|4.45l> holds with an additional 
factor , see (|4.73l) . that plays no role in the applications of this estimate. □ 


26 







































4.6. Local statistics of . In this subsection, we show that the gap statistics of the localized 
reference measure uit-l are universal, i.e., are given by the statistics of the Gaussian invariant ensemble 
up to negligible errors for large N. The precise universality statement for ujt-l is as follows. 

Theorem 4.10. There is a small universal constants e, %, a > 0, such that for any y £ Qt x (see (14.41 ) ), 
for any fixed, j and for any smooth compactly supported function O of n variables, we have 


E ^O^NqtM) (x i0 - * i0+j ))J = J = E G o((Vp # ) (xp o - H+i))J =1 ) 

+ 0(||0 / || 00 JV-«), (4.74) 

for N sufficiently large, for any io,i' 0 £ N n satisfying |zq — L\ < N x , \i' 0 — L'\ < N x with any 
L' £ [aN, (1 — a)N], and where g # := Q S c(Tl',sc) denotes the density of the semicircle law g sc at the 
location of the L'-th N-quantile of g sc . 

In short, Theorem 14.101 assures that the gap statistics of the localized measure ujt, in the bulk 
is determined by the Gaussian invariant ensemble in the limit of large N. This result follows from 
Theorem 4.1 of [28] and the properties of ujt, established in this section so far. 


Proof. Theorem 4.1 of [28], as stated, directly compares two local measures, but together with Propo¬ 
sition 5.3 in [28] it can also be stated as a direct universality result: if the conditions of Theorem 4.1 
of |28| hold for a local measure, then it has universal local gap statistics. 

Theorem 4.1 (see also remark after Lemma 4.5) in [2S] has two types of conditions. 

(1) Regularity of the external potential in the sense of Definition 4.4 of [28] . For the case at hand, 
the external potential V 1 defined in (14.23jl is regular if 

(vyw = ex, (x.) log (4kL) , (4.75) 

> -^y , x £ J z = [z_,z+], (4.76) 


hold, with some fixed constant c, where z 0 = {z+ + Z-)/ 2, d(x) = min{|a: — z+ 1, \x — Z-\} and d±(x) 
as in (I4.49[) . Here we used the notation Z- = zl-k- i, ~+ = Zl+k+i- 

In proving (14.471) with external points 7, we already established (14.751) . The convexity esti¬ 
mate (14.761) follows from the rigidity of 7: there is a constant c > 0 such that 


(V^)"(x) 




(2) The second input for Theorem 4.1 of [28] is 


E“ Tl a; ? ; — a,- < C- 


N 


i £ /, 


(4.77) 


where (a*) denote the /C equidistant points in J z . We have already established (14.771) in (14.681) . 

Based upon these two inputs, Theorem 4.1 of [2S] implies (14.741) . □ 


5. Universal gap statistics for small times 

In Section 2] we showed that the equilibrium measure wji of the dynamics (14.211) has universal gap 
statistics. In the present section, we compare the gaps of the two dynamics (14.211) and (14.201) . We 
proceed in three steps that are outlined in the Subsections 15.11 IST^l and [5751 In Subsection 15.41 we then 
complete the proof of Theorem 12.11 

As in Section H] we will fix a Y £ Q, or equivalently Y £ Q, but do not always indicate this choice 
in the notation. All estimates obtained will be uniform on Q, so that we can integrate out Y at the 
very end of Subsection 15.41 
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5.1. Step 1: Small scale regularization. First we introduce a small regularization in (14.201) starting 
from the time T\. This regularization is only needed for the critical case (3 = 1, where the level 
repulsion, c.f., Assumptions (3) of Theorem 12.11 is weakest. Level repulsion and the regularization 
introduced below will allow use to bound the kernel Bij defined below in (15.81) as E|i?, 7 1 < N\ see (15.141) . 
For (3 > 1, a similar bound may be obtained without any regularization. In the following we carry the 
regularization along since the case [3 = 1 is the hardest. 

This regularization procedure is the same as in Section 3.1 of da, but it is different from the regu¬ 
larization in the wji measure and in the x dynamics explained in part (2) of the proof of Proposition l4.9l 
(where the regularization parameter was called £*). Choose 


£jk 


e ■ l{j,ke i a ) 
-£ ■ l(j,k £ I a ) 


if j >k, 
if j <k, 


with £ := N~ 10Cl , 


for a large Ci > 1. (Note that by the above choice £* -C £.) 
Define the regularized version of (14.201) as 


(5.1) 


dxi(t) 


2 dBi(t)-v L dt+^J2 = 


/3N 


N Xi ^ ~ x ^ t ) + £ ii 


d t 


1 1 _ Xi(t ) + V k t ^ 

N Xi(t) — y k (t) + £j k 2 


i £ I , t £ [Ti,t 2 ] > (5-2) 


with initial condition x(Ti) = x(Ti), where the Brownian motions ( Bi) are the same as in (14.201) 
and (14.211) . Note that x may not preserve the ordering of the particles, but we will not need this 
property. 

Lemma 5.1. Define the event 


*== n 

t£[Tl,t2] 


< max | Xi(t) 
( i&I 


Xi(t) | < N 



(5.3) 


Under the conditions of Theorem, \2.1[ especially the level repulsion assumption (|2.19[) . there is a set 
Q* C Q with ¥(Q*) > 1 — N~ Cl such that P Y (S 1 ) > 1 — CN~ Cl holds for any Y £ Q*. In particular, 
the local statistics ofx.(t) and x(t) are asymptotically the same for any t £ [Ti,t 2 ]. 

Proof. Let 1Z be the rigidity set 

n := (|®i(t) -Ji(t)\ < Nt/N : t £ [T u t 2 ] ,*€/}. 

First we claim that P(JZ fl S 1 ) > 1 — N~ 2Cl . This estimate can be proved following the argument in 
Section 3.1 of da for vl = 0. Mutatis mutandis the same proof applies for vl 0. As an input, we 
need a level repulsion estimate of the form 

E1(R W+ lW -*, W+ Ell* - ]VW|: '° 8£| -’ Vte|r '- i21 ’ (5 ' 4) 

that follows from (12.191) . Therefore, by conditioning, there is an event Q such that P (Q ) > 1 — N~ Cl 
and _ 

p Y (KDH 1 )>i-/r Ci , VYeif. 

Since P (Q) > 1 — N~ D for any D > 0, see (14.7|l . without loss of generality we can assume that Q C §. 
Note that for Y £ Q we have P Y (7?.) > 1 — N~ D , for any large D > 0; c.f., (14.51) . Choosing D larger 
than Ci, completes the proof. □ 
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(5.5) 


5.2. Step 2: Holder regularity. To compare the gaps of x and x, we introduce 

v = v(i) := e (‘- T 0/ 2 (x(t) - x(t)), t > T[. 

Subtracting (15.21) from (14.211) and dropping the t argument for brevity, we have 


d Vi 1 ^ 

~dt “ ~~ ^ 






____ ____ i_____—e^ ( — Ti) 

N (Xi -Xj +£ij)(xi -Xj) ~‘N ££(xi-y k + £ik){xi-%) 2 

j¥=i 


+ J_ (xj - Xj) - (xj - Xj) + £jj + J_ {Vk ~ 7fc) + (%i - x t ) + £ ik 
N fti (Xi-Xj +e ij )(x i -x j ) N (Xi-y k +£ ik ){xi- 7 *,) 
j¥* 

i £ I. We rewrite (15.61) in the form 


(5-6) 


^ = -(^v) j + t; (1) + t; (2) , 


{SSv)i := ^ Bij(vi - Vj) + WiVi , 

ie-f 

with time-dependent (symmetric) coefficient^, € /, 

1 := — V —_-_ - 

N (t.a — II, 4- Pa 


. _ _ 

N (Xi Xj £%j) {X-i Xj ) 


N £gi ( x i ~ Vk + Sik)(xi - 7fc) 


and with the “forcing terms” 

P (l) J_ (Xj — Xj) - (Xj — Xj) + £ij _ 1 (t-T')/2 


F U) , = \Xi-Xj) , , ~ - _ _ e v 

N ( Xi — Xj + £ij)(xi — Xj) 2 


VL(t - Ti), 


F ( 2 ) , = IW _^ 




l 




24 - 7* 


N i^l( Xi ~ y k + £ ik){xi -7*0 ^ 0*7 - Vk + £** 0(27 - 7*0 


(5.7) 

(5.8) 

(5.9) 
(5.10) 


(Since 27 , 27,27 and y k depend on time, we have Bij = Bij(t), Wi = W»(i), etc.) 

We first study in Subsection 15.2.11 the “free dynamics” generated by £8 and then deal with the 
forcing terms (7*4), ( F.f 2 ' 1 ) with a perturbative argument in Subsection 15.31 

5.2.1. Holder regularity of the free dynamics. Let v solve without the forcing terms, i.e., 


du, ; 


-P- = ~(®v)i = - Bij(vj - Vj) - WiVi , t>T[, 


(5.11) 


with initial condition v(T[) = v(T{). 

We will need a certain upper bound on the coefficients Bij in a space-time averaged sense. Let 
T := \T{,T"\. Mimicking Definition 9.7 in [28], we say that the Equation (15.111) is regular at a 
space-time point (Z, 9) £ I xT with exponent p > 0 if 


sup 

ter 


sup W~ 

\<M<K JV 


l*- 0 | 


f Jf Y, J2 \Bij{s)\ds < N 1+ P. 

Jt ieI:\i-Z\<M j£l:\j-Z\<M 


(5.12) 


Furthermore, we say that the equation is strongly regular at a space-time point (Z, 9) £ I x T with 
exponent p > 0 if it is regular at all points { Z} x {0 + D}, where the set Q is defined as 

fi | ■ 2 -m (l + 2~ k ) : 0 < m, k < Clog A'} . 

From Theorem 10.1 of [28] and Lemma PTTI we obtain the following result. 

TSometimes we write /i, , instead of Bij to clarify the notation. 
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Lemma 5.2. Let c\ ~ 1/100 and choose T" = T[ + K Cl /N. Then, there is an event z? and constants 
C and p = p(i ;) > 0 such that on the event S 2 the Equation (15.1111 is strongly regular at ( L,T"), and 

l(H 1 nH 2 )|tT i+ 1 (T")-^(T")|<C'iV- 1 +^-‘’/ 4 , \i — L\<C, (5.13) 


where q > 0 is a universal constant. Moreover, we have the estimate P Y (S 1 n^ 2 ) > 1 — N p / 8 , for N 
sufficiently large. 

Proof. We apply the Holder regularity result, Theorem 10.1 of [28] . to the evolution equation (15.111) . 
Thanks to the regularization introduced in Step 1, c.f, Subsection 15.11 we have, for any i,j G I and 
t G [T[, t 2 ], that 


E Bij < N ( maxE—-— 

3 - V i N \Xi - Xi- 1 + e|]p 


1 /p / i \ 1 /q 

( maxE-—-—-7, 

\i&I [N\Xi - Xi-^ / 


~ N (Ne^lNlxi-Xi.x+e |] 2 

<N(£-*N s+ Z+*\\oge\) 1 /p C qW 
< CN 1+5+2/i , 


i/p 


maxE-—-—-n - 

V [N\xi - Xi-xW* 


1/9 


(5.14) 


where we first applied Holder’s inequality with conjugate exponents p, q , with p = 2 + </> > 0, then 

used (15.41) and (14.461) . and finally chose <p sufficiently small depending on C i in (15.11) . 

Notice that to guarantee regularity in the sense of (15.121) (modulo a constant factor), instead of 
taking suprema over all s G T, M G Jl,/CJ, it suffices to take suprema over a dyadic sequence of times 
Sfc = 6 ± 2~ k and parameters Mi = 2 l , k,l G J1, Clog/VJ, since space-time averages on comparable 
scales are comparable. Using (I5.14L setting p := 6 + 3£ and applying Markov inequality, for any 
fixed values of s and M , the space-time average in (15.121) is bounded by N 1+p with probability at 
least 1 — N~ p / 2 . Taking the union bound for not more than C(logiV) 2 times, we can guarantee 
regularity at any space-time point with probability at least 1 — N ~ p / 3 . Since the definition of strong 
regularity requires regularity at not more than C(\ogN) 2 space-time points, an additional union bound 
guarantees strong regularity at any fixed space-time point with probability at least 1 — N~ p / 4 . Defining 

S 2 := {Equation (15.111) is strongly regular at (L, T[') } , 


this proves that P Y (S 2 ) > 1 — N ~ p / 4 and verifies condition (Cl) in Theorem 10.1 of [28| on S 2 . 

The other condition (C2)^ in Theorem 10.1 of [2S] concerns large distance estimates of Bij. More 
precisely, condition (C2)^ requires that 

TV 1-5 

te[T[,Tf], (5.15) 

I* J I 

for any i,j with \L — i\ < K/C , |L — jj < K/C, and that 


iVl(min{|L - i|, |L - j | } > K/C) < < _ CN 


C\i~j\ 


I* - 3\ 


t&[T[,Tf], 


(5.16) 


for any \i — j\ > CN with some constants C,C > 10. Further, W, is required to satisfy 


CN N CN 1+ t 
<W l(s) <^- 


t G [T{,T"] , 


(5.17) 


where A, := min{L + K + 1 — i, L — K — 1 — i}. Using the rigidity estimates for x, x of Lemma 14.91 it 
is easy to check that, for any £ > 0, there is an event S 2 and constants c, with P Y (S 2 ) > 1 — e _cJvS , 
such that (15.171) . (15.161) and (|5. 151) hold. Set S 2 = S 2 D S 2 . Then for all sufficiently small £ > 0, we 
have P(S 2 ) > 1 - CN~ p / 4 . 
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Let HAIloo := sup ig7 \Ai\, A £ C N . Then the conclusion of Theorem 10.1 of [28] for the equa¬ 
tion (15.111) is that 

l(H 2 )|u i+1 (T")-^(Ti")| ^CA-^llv^JIloo, I i~L\<C, (5.18) 

where q > 0 is a universal exponent and where T" = T-[ + K Cl /N. More precisely, (15.181) follows 
from (10.6) of [28] after rescaling space and time by setting the constant a equal to 1/4 (this a is 
different from the a used in the present paper). 

Next, recalling from (15.111) that Vi(T[) = Vi(T{) and that Vi(T {) = Xi(T{) — we get 

MTDWoo < PW) - X^OIU + ||x(T() - STOIU < CN~ 5C 0 + CN~ 1+ * , (5.19) 

on 3 1 fl S 2 , where we used Lemma l5Tl and that the processes (x(f)), (x(t)), are both rigid in the sense 
of Definition 14.81 for t £ \T\, t-f\- Thus, combining (15.181) with (15.191) we get 

1 (S 1 D H 2 )|ujj_|_ 1 (T") — Vi{T")\ < CN~ 1+ ^K~ q / 4 , I i-L\<C, (5.20) 

where the event E^flS 2 satisfies P y (S 1 DE! 2 ) > 1 —lV _p / 8 , p = p(ff) > 0, for sufficiently small £ > 0. □ 


5.3. Step 3: Removing the forcing terms. Having established the Holder regularity of the free 
dynamics of (15.111) . we now deal with the “full dynamics” (15.61) . The main result is as follows. 

Proposition 5.3. Let ci ~ 1/100 and choose T" = T[ + K Cl /N. Then there is an event 3 and 
constants C and C 2 , C 3 > 0 such that, for any y £ Q, 





N ~ c 2 

N 


(5.21) 


for N sufficiently large, where jl := \L — K/A,L + K/ 4]. Moreover the event 3 is such that 3 C 
and satisfies P(S D Q) > 1 — N~° 3 , for N sufficiently large. 

The proof of Proposition 15.31 is given in the following subsections. 

5.3.1. Removing the forcing terms F^ . Subtracting (15.111) from (15.71) we have 

d{Vi - Vi) =-[mv-v)\.+ F i , i £ I, t£[T{,T{'}. 

Eventually, we will choose i £ jl := \L — K/4, L + A//4J, yet here we can take i £ I. 

Let ^ag(t,s) denote the time-dependent propagator of the equation (15.111) from time s to t, with 
s < t. From Duhamel formula we have 


i{t) - Vi(t) = [ {Wgg(t, s)F(s)) i ds , 
Jt> 


i £ I , t>T[. 


Note that is a contraction in the sup norm by maximum principle (recall that W) > 0). Thus 


i(t)-Vi(t)\< f max|T) ( 1 ) (s)|ds + f (tygg{t, s) F (2) (s)) 


ds, 


t>T[. (5.22) 


Fixing t = Tf, using Lemma T5. II and the choice of e in (15.11) . we estimate 


i(5 1 n 5 2 ) 


IT' 


max ^(s)|ds < 1 ( 
iel 


-—'2 


)CN 


—5Ci 


IT' 


EE 

i£l j£l 


Bijis^ds + CvLiT? -Tr) 2 , (5.23) 
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where we estimated the maximum by the sum. Thus recalling (15.1211 and using (15.231) . (15.221) we get 


1(S X DS 2 ) [ 1 max |Fj ( 1 ) (s)|ds < CN~ 5Cl K(T" - T[)N 1+P + Cv L (T" - Ti ) 2 
Jt[ ie/ 

K i+ Ci 

< CK 1+Cl N~ 5Cl+p + c . 

N z 

Since C\ > 1, we conclude that the effect due to F^> is below the precision we are interested in, i.e., 
there is c > 0 such that, for any i £ /, 


i(s x ns 2 )|u i (T") 


< civ - 1 " 0 + i(s 1 n s 2 ) 



(^(T", a)FW(s)). 


ds. 


(5.24) 


5.3.2. Removing the forcing terms F ^ 2,ln \To estimate the influences of the forcing terms 
write 


(Fi 


(2)x 


p( 2 ) _ p(2.in) _|_ p(2,out) 

r i 1 i ' r i i 


with 

^ (2 ' in) := if } 1 (i £ \l) , if’° Ut) := if > 1 (i £ / , » 0 ±/) , 

where we introduced |i := [L — TT/2, L + if/ 2 ]. 

To control the inside part if ,in \ we use the following lemma. Recall the definition of the event Q 
in (14.51) and the definitions of the intervals of consecutive integers Iq and I a in (14.ip . 

Lemma 5.4. Let K satisfy (|2.27l) and fix Y € Q . Then we have the following estimates. 

( 1 ) For all k £ I<r\I, we have 

\y k (t)^y k (T 1 )\<C^-+CN s (t-T 1 ) l -^^-+CN s (t-T 1 )\ i£[Ti,i 2 ]. (5.25) 


(2) For all k £ Jo\/, we have 




t £ [Ti, tf\ ■ 


(5.26) 


We complement Lemma 15141 with the estimate 

\y k (t)-j k \<CN^Vi, t>T i, k€l c a , (5.27) 

on Q, as follows immediately from the definition of Q in (14.51) : c./.. Assumption (4) of Theorem 12. II 
Proof. To prove (15.251) . we estimate 


\Vk(t) - Vk( T i)l < I Vk(t) - 7fe(*)l + IlfcO) - 7 fc(Ti)| + \y k {Ti) - 7fe(Ti)| 

on the event Q, where we used the rigidity bound in (14.21) for k £ I a . Next, we write 

7kCO _ lk{Ti) = Tk{t) ~ v L {t - Ti) - 7 fe(Ti) = f %(s) ds - v L {t - . 

JT ! 


(5.28) 
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Then by Lemma 14.31 we have 


[ ik(s)ds= / 7 i (s)ds + 0(JV- 1 + 5 (t-ri)|fc-L|) 

JTx Jt-l 

= 7 iCTi)(t - Ti) + 0(A^- 1+5 (i - Tx) 2 ) + 0(iV 1+5 (f - T x )\k - L\). 

Thus recalling that Vl = 7 l(Ti) by definition, we conclude that 

| 7 fc(i) - 7 fc(ri)| < CN~ 1+ \t - Ti ) 2 + CW" 1+4 (t - T\)\k - L \. (5.29) 


Together with (15.281) this implies (15.2511 . 

To bound the left side of (15.261) . we split 


I2/fcO) - 7fe| < I Vk(t) ~ 7k(*)l + l7fe(*) - 7fcC r i)| + |7fe( r i) “ 7fc| 


(5.30) 


Then, using the rigidity from the definition of Q\ in (14.21) . the first term on the right side can be 
bounded by CN^/N. The second term on the right side is controlled by (15.291) . To bound the third 
term on the right side we apply Corollary 14.71 to find 

l7fe(ri) ^ 7kI < C^- + c N * hk K lLl + CN S \j k - 7 l| 2 . 

Recalling that |7 k — jl\ < CK 5 /N, for k G Iq\I , and using that K satisfies (12.271) . we get (15.261) . □ 

We now bound the term F^ 2,1D \ Abbreviate 

1 


Bik — 


N(xi - y k + e ik ) (Xi - 7 k ) ’ 


* g i , k g r. 


(5.31) 


Recall the bound on l(2 1 )|a;i — Xi | from Lemma [5.II and the definition of (e^) in (15.11) . Using 
Lemma 15.41 and recalling that t — T-\ < CK 2 /N, we obtain 


1 (s 1 )!^ 


(2,ii 


E 

kei 0 \i 


N5G 1 


Nt ( Nt\k-L 

~N + 


N* 

~N 


■c E 

\/ 0 

c v ( NtJ r= 


■ N' 


K N 

S K 2 \L — k\ 


N N 
l-Bzfc(s)|, 


Bik(s) | 


|-Bjfc(s)| 




s G Pij T"} , (5.32) 


where we also used (15.271) together with T" — Ti< CK 2 /N to get the last term on the right of (15.321) . 

To perform the sums over k in the first two terms on the right, we recall (15.311) and we note that 
there are two constants c, d > 0 , such that 


\Vk ~ x i\ > C \L ~ k\/N , 


17k - Xi\ > c'|L - k\/N , 


k G I a \I , 


(5.33) 


where we used the rigidity estimate for y (embodied in Q\ see (|4.5I) L the rigidity estimate for 7 obtained 
in Corollary 14.71 and the rigidity estimate for x, x obtained in Proposition 14.91 as well as the choice 
i G \l = \L — K/ 2, L + K/ 2J. The summation over k qL I a in the third term is estimated using that 
| y k — Xi\ I 7 k — Xi | > c"(cr) > 0, k £ I a - Hence, after summing up the right side of (15.321) . we get 


l^lif’ in) ( S )| < 


1 


< C- 


N^-ir 

iV 2 ? 


Ni N 2 t 7V« 

K + K + NK 3 


C 


N^K 

Vn 


K 




s G [T [, f 2 ] 


(5.34) 
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where we used 5Ci — 1 > 1 and that K satisfies (12.271) . Thus by (I5.34|) we have 



F i (2,in) (s)| ds < (T" 



K Cl 7V 2 « 

~W~K~ 


for all i £ \l 7 so this error is below the precision we are interested in: For some c > 0 we have that 


N~ c 

lts 1 )!^') -^(T")| < C7— + MS 1 ) 



* e 2 J - 


(5.35) 


The outside part F < - 2,out ' ) is treated with a finite speed of propagation estimate. 

5.3.3. Removing the forcing term F^ 2,out \We first recall the finite of propagation estimate for the 
propagator Abbreviate for simplicity ^(t,s) = c flfgg(t,s) and denote its kernel by ij{t,s ), 

i.j £ I. By Lemma 9.6 of [33] there is C such that 

K x l 2 \lN(t - si + 1 

s)| < C- V ,. E , i,j€l, t > s >T\. (5.36) 

on S 2 . We refer to (15.361) as a finite speed of propagation estimate. 

Next recall that we want to control 

max U ij (*> s)Fj 2 ’ out - > = max E E E (t, s)F ) (2) (s), 

j&i 

where T[< s <t < T" , and where we have introduced 

jA 2 )( s ) % ~ Fi + _|_ ~ \ / n _ 

’ V (Xj - Hfc + Ejfe)(ajj - 7fc) (Xj ~Vk + £jk)(xj - Ik) J 

With some large C, we next split the summations over k and j as 


s)F^ out) = 1 E E 1 (lj - fc l > CNt)Uij{t, s)F$\s) 

jei j e i\^ikei a 

+ ^ E E 1 ^- fc i <c '^)^( t ’ s ) ii, ifc ) ( s )> iG i 7 - ( 5 - 37 ) 

je/ \i/fce/<= 


1 

iV 


We start with bounding the first term on the right side of (15.371) . On the event S 1 , we can bound 

1 Nt N*\k-L\' 

HI J - O ^ ^ JV 7 ( jy 

fce/o\/ 


fce/ c kei 0 \i v 2 


Fjfe(s)| 


+ c E Mli-*l>c»E)(^ + w^!LA! 


kei„\i 0 


I-®a( s )I 


+ c E i(li - ^1 > CN ^^§- \ B M\ ■ 

k£I(j 


(5.38) 


here we used (15.361) and the Lemmas 15.11 and 15.41 We further used that s < t, t — Tf < CK 2 /N by 
assumption. Thus, summing over fc, we get 


i(s x )^E 


fce/° 


< c 


AT 2 ? 


|j — L + K + iV«| ’ 
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where we used the estimates in (15.331) . See (15.321) and (15.341) for a similar estimate. Returning to (15.371) 
we see that the first term on the right side is bounded as 


N £ £*( \j-k\>CN^\U ij (t,s)F^(s)\<C E 

jei\UWi 


N 2 tK 1 !' 2 ^N(t - s) + i 


^ \j-L + K/4\\j-L + K + Nt\ 

j : K/2< j 1, <K 1 UJ 1 

< CN 2i K - 1/2+C1 / 2 , (5.39) 

on S 1 nS 2 , where we used that t — s < I\ Cl /N. 

It remains to control the second term in (15.371) . Similar to (15.381) . we have on H 1 , for j £ I \ 4/, 
that 

i^|i(|,-t|<c»£)4 2 »w|<c V 1(|3-*I<CM)( * +^ + ^i) |S*WI 

fce/ c kei 0 \i v ' 

E H\j-k\<CNt)\B jk (s)\. 


kei 0 \i 


-£ 

N ^ 

kei c 

< 


A ± ’ 
rT ' 1 


We thus have, for i £ jl, 

ds 1 (|j - k\ < CNt) Uij (T" a) F$ (s) 

c w C ds E E Kb- - k < |S j4 wi 


< C 


k£l c j£l 

N 2 ^ A ' 1 / 2 / ’ T 


L+Jf 


jtttI d. y My+y>i 

V 17 T 1 y = L +, ff -| CW«J 

^2^1/2 _ /-R" ^-^+L CJV 'l 


-C- 


\/NK 


, Tl " OiV | 

I ds E I^e-i( s )I 

■' T ' i=L-K 


(5.40) 


where we used that \Bjk\ < I? & > j, respectively |IL,- fc | < k < j, for all k £ I c , j € J. 

(Here and below we use the convention that, for j £ /, Bj L ±^ K+ 1 ) = R^z.jyif+i)! respectively Bj t k = 0 
if |A;| > L + I\ + 1.) To bound the two terms on the right side of (15.401) . we use that the evolution 
equation (15.111) is “regular” at the space-time points (L + K, T") and (L — K , T”). 

Lemma 5.5. There is an event S 3 such that evolution equation (15.111) is regular at the space-time 
points ( L + K , T") and ( L — K, T") in the sense that 


l(^ 3 )sup sup — 

sET1<M<K. bv 1 + |S 


T"| 



E |-Bj,»±i(s)|ds < N 1+p . 

iel ■■ \i-L±K\<M 


(5-41) 


Moreover, we have the estimate P Y (^ 3 n S 2 D S 3 ) > 1 — N p / 10 . 

Proof. We can follow almost verbatim the first part of the proof of Lemma [ 5 ~ 2 l Using (15.41) and (14.461) . 
we can bound E|-Bi,j±l| as in (15.141) . Then dyadic decompositions around the space-time points (L ± 
AT, T") combined with applications of Markov inequality yield the claim. □ 


Next, returning to the estimate in (15.401) . we conclude from Lemma [53] that the first term on the 
right can be bounded as 


i(s 3 ns 3 ) 


N^K 1 / 2 

Vnk 


y/T" - T l 


L+K 

d S e 

j=L+K-YCNi\ 


N 3 i 

B jJ+ i(s )| < Cj=(T[' - T[) 3 / 2 N 1+P 


< CN 3 ^K 3ci ^ 2 N~ 1+p . 


35 













































Using the same argument to bound the second term on the right side of (15.401) . we conclude that 


nE3 )jvS 


fce/° 


ds 1 (| j -k\< CN^)Uij(T", s ) F%\s) 


jk 


Kr3^+p K 3 Cl /2 i 

---, i £ —I. (5.42) 


Summarizing the estimates above, we can now state the proof of Proposition 15.3 


Proof of Provosition \5.3\ Let E := E 1 HE 2 DS 3 . Note that P(S) > 1 — N~ C2 , for any 0 < C 2 < p (with 
p = 6 + 3f). Adding up the estimates (15.421) . (15.351) and (15.241) . and recalling that K satisfies (12.271) and 
that ci ~ 1/100, we conclude that there is C 3 = C 3 (£) > 0 such that (15.211) holds, for £ > 0 sufficiently 
small and N large enough. □ 


5.4. Conclusion of the Proof of Theorem 12.11 Recall from (15.51) the definition of (i>j(f). Combin¬ 
ing (15.211) and (15. 13|) we obtain 

l(Z)\v i+1 (t 2 ) - Vi(t 2 )\ < N- 1 - 0 *, \i — L\ < C, (5.43) 

with some small C 4 > 0. Moreover, we have P(SD^) > 1 — N~ C5 , for some C 5 > 0. This exactly proves 
the following result. 

Lemma 5.6. The gap statistics o/x(T") andx(T") for indices near L coincide in the limit of large N. 

Combining this with Lemma 15.11 we need only to understand the local statistics of x. But by 
Lemma 14.51 this is the same as the gap statistics of the local equilibrium measure uit 0 ■ The latter one 
is universal as we showed in Subsection l4.61 To conclude the proof of Theorem 12.11 we note that we can 
integrate over Q, as follows from Lemma 14. II and the assumption that the observable O is compactly 
supported. Finally, choosing Tf > t\ such that T" = T, we obtain (12.211) . This completes the proof of 
Theorem 12.11 


A. Semicircular flow 

In this appendix we study the semicircular flow in more detail. In Subsection lA.il we prove Lemma 14.21 
and Lemma 14.31 In Subsection IA.2I we discuss the Assumption (1) of Theorem 12.II in more detail by 
arguing that it is satisfied for a large number of random matrix models. 

A.l. Classical flow of the density.Recall from (12.61) that rn t satisfies 

m t (z) = [ — m - g(2/)dy — -, lmm t (z)>0, lmz>0, (A.l) 

Jr e~ t/2 y - (1 - e- t )m t (z) - 2 

for all t > 0, and that rnt determines a density g t via the Stieltjes inversion formula, i.e., gt(x) = 
L lim^x^o Im mt(x + i?;), i£l. We call the map 1 1 —» gt the semicircular flow started from g. 

It was shown in { 8 ] that the density g t is a real analytic function inside support for fixed t > 0. 
Yet, without any further assumptions on g , estimates on the derivatives of g t deteriorate for small f, 
since the Equation (IA.1I) may lose its stability properties (i.e., the denominator on the right side can 
become singular). This can, for example, be remedied by imposing the conditions in Assumption (1) 
of Theorem [241 Consider for E > 0 the domain 

:= {z = x + ip £ C : x £ [E — E, E + E], 77 > 0} . 

Denote by mo the Stieltjes transform of g. In accordance with the Assumption (1) of Theorem 12.11 
we assume here that mo extends to a continuous function on T>s and that there is a small 6 > 0 and 
a constant C such that 

sup |too( 2 )| < C, sup \d%mo(z)\ < C(N s ) n , n=l,2. (A.2) 

zGX> z£X> 
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Lemma A.l. Consider the semicircular flow gt started from g. Let g satisfy Assumption (IA. 2 I) with 
exponent 8 > 0 and £ > 0 . Then there is C' > 0 , such that 7774(2) is uniformly bounded on T>y, for all 
0 < t < CN~ 2S . Moreover, there are constants C,C', depending only on g, such that 

sup | m t {z) — 7710(2)1 < CtN 6 , (A. 3 ) 

2 


for all 0 < t < C'N 2S . Further, there are constants C,C’, depending only on g , such that we have 
the bounds 


sup \dfmt(z)\ < C(N s ) n , n = l, 2 , 

Proof. Set 04 := 1 — e _t . Starting from (IA. 1 I) we obtain, for t > 0, 


0 < t < C’N - 25 . 


(A. 4 ) 


\ m t( z )\ 2 < 


Q{y)&y 


|e t / 2 y - z - a t m t {z )\ 2 J ij + a t lm.m t {z) ’ 


11717774(2) 


G C+ 


where we first used Schwarz inequality for the probability measure g to get the second line. Then we 
used once more (jA.ll) . We thus obtain the rough a priori bound 


\m t (z)\ < a t 1/2 , 


t > 0 , 


for zeC + Ul. Next, we introduce 


m t {z) ■■= / 
Jr 


g(e t / 2 v)e t / 2 dv 


v — z 


(A. 5 ) 


(A. 6 ) 


Note that nit is uniformly bounded on, say, T>s/2 for, say, t < 1 . This may be seen by writing 
fh t (z) = e t ^ 2 mo(e t ^ 2 z). Then we can write 

m t (z) = m t {z + a t m t (z)). 

Thus, using the estimates on mo in (IA. 2 I) . we have 

| mt(z) - 777 t( 2 )| = (7774(2 + a t m t (z )) - m t (2))| < CN S a t \m t (z)\ < CN s a ] /2 , (A. 7 ) 

0 < t < 1 , on T>s/2, where we used the a priori bound (IA. 5 I) . It follows that 


T — 2 6 


(7774(2)! < \m t {z)\+CN 5 a ] /2 <C, 0 <t< CN~ 

B7.1t then reasoning once more as in (IA. 7 I) . we must have 

|m t (2) - 7774(2)1 < 17774(2 + 0-47774(2)) - 7774(2))! < Ca t N s \m t {z)\ < Ca t N s , 


(A. 8 ) 


-25 


0 < t < CN , for all 2 G T>s/2- We hence obtain that (777.4(2)! < C on Vy, and 0 < t < CN 
Next, we observe that 

1 7774 ( 2 ) - too( 2)| = |e t/2 ?77 0 (e- t/2 2) - mo(2)| < C{e t/2 - 1) + CN s (e t/2 - 1). (A.9) 

Combining (IA.8I) and (IA.9I) . HA. 31) follows since t < C'N~ 2S by assumption. 

To deal with the derivatives of mflz), we first note that we for 2 G T>y/ 2 , we have e^^+o^rm^)) G 
T>y, for t < 1. Thus we can bound, for 2 G T>y/2 an d t < N ~ 2S , 




g(v)dv 


= a t \(d z m)(z + 0-4777.4(2))! < CN s a t < CN 5 , 


(e t / 2 v — 2 — at.mflz )) 2 
where we used the definition of 777(2) in (IA. 6 I) and the assumptions in (IA. 2 I) . 


(A. 10 ) 
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Next, differentiating (IA.1I) with respect to z, we obtain 


g(v)dv 

( e~ t / 2 v — z — atmt(z)) 2 


{d z m t {z)) 



(e _t / 2 v 


g(v)dv 

- z - a t m t {z)) 



Hence, using twice (IA.10I) . we get \d z m.t(z)\ < CN S , for z £ Vy ,/2 and 0 < t < CN~ 2S . Repeating the 
argument, we see that there is another constant C such that \d 2 rrit(z)\ < CN 2S , 0 < t < N ~ 2S , for all 
z £ 2? s / 2 . This proves (IA.4I) . □ 


A.1.1. Quantiles. From (|A.4[) . we see that the derivatives of mt.{z) are bounded inside X>s /2 for 
t < Cl. Without further assumptions on g we have little control on mt(z) outside T>y./ 2 - Yet, we can 
circumvent this problem, by introducing a regularization of mt(z), respectively the measure gt, as 
follows. Throughout the rest of this appendix, let 77 * > 0 satisfy 

V*N 5 « ^ • (A.11) 

Recall the definition of the Poisson kernel P. in ( 12 . 21 ) . We then set 

qT(x) := (P^ * g t ){x) = -Imm t (a; + b?*). (A.12) 

7 r 

We claim that 

eT(y)&y , . . A 

-= ttjRz + 177 *), 2 €L t . (A.13) 

y - z 

Indeed, using (12.31) and (IA.12I) . we have 

-Im [ Qt = {Pr, * gT)( E ) = ( p v+r ). * et){E) = -Im m t (E + irj + i? 7 »), 

where we used P v * P^ = P^+tj, ■ Since z = E + ir] and since the Stieltjes transform is analytic 
in the upper half plane, we get (IA.13I) . In the following we write ml* (z) := nit(z + i? 7 *). Note 
that gl* is a probability measure. It follows from basic properties of the Poisson kernel that gl* 
converges uniformly on compact sets to g t as 77 * \ 0. Using (IA.2I) it is then easy to check that 
\Qo*( x ) — PofY)! < CN s r]x -C N -1 , for all x £ [U* — E,U* + £]. Since the semicircular flow preserve 
regularity (for short times see Lemma [A.II) . we also get \gl*{x) — Qt(x)\ < CN 5 i 7 * -C N~ 1 , for all 

x e [£*- e/ 2 , u* + s/ 2 ], o < t < n~ 2S . 

As a consequence of the regularization in (IA.12I) . gl* is smooth with bounded derivatives (in terms 
of inverse powers of 77 *) that all lie in L P (R), 1 < p < 00 . Consequently, the following basic properties 
of the Hilbert transform can be justified easily (see, e.g., |47j): For n £ N, 

d n (Tq?) = {-l) n (T(d n gl*)) , (A.14) 

(here d n denotes the n-th spatial derivative). Further, we have T (Tg/*) = — g /*. 

Next, we define the continuous quantile, 7,„,* (t), w £ [0, AT] , of the measure gl* by 

nZ* (*) 

/ g?(y)dy=%, e V o*=Q v *, (A.15) 

J —OO ^ 

c.f., (12.91) . Note that 7 lf(t) is defined for any w by (IA. 151) . since the measure gl* is supported on 
the whole real axis. The measure g t (without regularization) may be supported on several disjoint 
intervals. This leads to some ambiguities in the definition of the quantiles, c.f., (12.91) for one way of 
resolving the ambiguity. Using the regularized density gl* is another way of avoiding this ambiguity. 
Nonetheless, we emphasize that the ? 7 *-regularization is simply a technical tool: for every practical 
purpose we have 77 * = 0 and the reader may forget about it in the subsequent arguments. 
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Corollary A. 2. Under the assumption of Lemma \A.1[ the following holds. For 0 < t < C'N 2S , 
with C' sufficiently small, we have the estimates 

| gf ( E) - g" ( E )| < CN s t , | {TqT ) (E) - (Tg r,r ){E) | < CN*t , (A. 16) 

for all fiEi. Further, for w such that J^*(0) £ [E* — £,!?* + £], we have the estimate 

l7S*(i)-7S*(0)|<CA- 5/2 , 0 <t<N~ 2S , (A.17) 


for some C. In particular, if 7 ^,*( 0 ) £ [E* — S/4,1?* + S/4], then 7 %f(t) £ [1?* — S/2,1?* + S/2], for 
all 0 < t < C'N~ 25 . 

Proof. The estimates in (I A. 161) follow from (1A.3I) by noticing that y/*(l?) = n~ 1 Immt(E + iry*), 
respectively (Tg/*) (1?) = Rem t (l? + i?j*). To establish (IA.17I) . we note that, by definition, 



eT ( V ) d y 


Thus, using (1A.16I) . we get 


r'Yw* (°) w 

/ Q’ n *(y)dy = -. 

J — OO ’ 


'(*) 


(y) d y = 



g r> * (y) d y + 0{y/Wt) , 


(A.18) 


where we also used that g has finite second moment. By our assumption on w we must have 
0 v *(7w(O)) — c > f° r some c > 0. We thus get from (1A. 181) and (IA.2I) that 



pZ* (*) 

pZ" ( 0 ) 

e’ , *(72*(o)M*W-72*(o)|<c 

/ Q rh (y)dy - 

/ Q v * {y)dy 


J — OO 

J — OO 


< C(N s t) 1/2 . 


Thus, for 0 < t < CN~ 25 , we get (t) - 'yff (0)| < CN~ 5 / 2 . □ 

The estimate (1A. 171) will serve as a priori bound below. To get precise estimates, we derive next 
the equation of motion of 7 ,/,* (f) under the semicircular flow t —> gf*. 

Lemma A.3. For t > 0, we have for all w £ [0, N], 


and 


d 7£*(*)_ r,,w 7 V* ( T ^*) (iwit)) 

dt MlwW) 2 2 gf" (jw(t)) 


^7w(t) = 1 

d w Ng^ilwit)) ' 


In particular, if gf* ( 7 ^* ( t )) > c, for some fixed c > 0, we have the uniform estimates 


d lw (t) 
dt 


= -M*)(7*•(*))- 


7w (t) 


0 (??*), 


d w 


= 0(N~ 1 ). 


(A.19) 


(A.20) 


(A.21) 


Remark A. 4. Lemma [A. 3 1 directly controls 7 %f(t) in the bulk. With some more effort the third term 
on the right of (IA. 191) can by controlled at the edges. For example, assuming that gt vanishes as a 
square root at, say, its lowest endpoint, it can also be shown that gf* ( 7 JJ* (f)) > y/rff, for small w. 
Thus the “error” term in (I A. 191) is of order yfrff, as ??* \ 0, at the lowest edge of the density gf*. 

Proof. We recall that mt(z), z £ C + , defined in (IA.1I) . satisfies the following complex Burgers’ equa¬ 
tion [48] (see also [55]) 


d m t {z) 
dt 


l_d_ 

2 dz 


( 'm t (z) ( m t (z) + z)^j 


z eC + ,t> 0 . 


(A.22) 
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Indeed differentiating (IA.1I) with respect to time we obtain (IA.22I) after a series of elementary manip¬ 
ulations. We use (IA.22I) with (z) = rrit(z + i/y*) replacing mt(z) in the following. 

To deal with the right side of (IA.22I) . we note that 


'(*) 2 = 


qT (aM-'e qT ( y) d v 

x — z y — z 


= 2 


r g ?* (apds eT ( v) d V 
Jr 2 x-z y-x 


= 2 


(T QT){x)eT (z)da; 


(A.23) 


z £ C + . Plugging (1A.23I) into (IA.22I) . we get 

dw?* (z) _ l_d_ f f 2 (T eT){x)eT (x)dx + f xg v t - (x)dx 


dt 


2 d z 


x — z 


x — z 


f i t 1*qT( x ) 

Jr x-z 


z £ C + , 


where we used that J R g^* (x)dx = 1. Differentiating the right side with respect to z, we get 

zeC+. 


d mf (z) = _1 I" (2 (T g'l”) (x) + x + iy*) eT (x)dx 
d t 2 Jo 


(x — zY 

Integrating by parts in x, taking the imaginary part and the limit ?/ 0 (with ?y* > 0), we obtain 


(>T{E) = \ [(2 (T 0 r) (E) + E) erm' + £ [«*) (E)}' 


(A.24) 


where we use the notation a = dta and a' = Secl, for any function a = a(t,E). On the other hand, 
differentiating the defining equation (IA.15I) of y2* (t) with respect to t. we get 


7 Z(t) = - 




e?* «*(*)) J-c 


et ( y) d v ■ 


(A.25) 


Hence, combining (IA.25I) and (IA.24I) we get (1A.19I) . 

Finally, to establish the first estimate in (IA.21I) . we note that (T g^“) is uniformly bounded by 
Lemrna lA.il Together with the assumption Pt* (lw* ft)) > c > 0, (IA.21I) follows. To prove (IA.20I) and 
the second estimate in (IA.21I) it suffices to differentiate (IA.15I) with respect to w. □ 


Remark A.5. For g £ A4(R), let 

Ent[e] := [ \x 2 q(x) dx - [ log \x - y\dg(x) dg(y) , (A.26) 

Jr z Jr 

Voiculescu’s free entropy. Then the limiting equation of (IA.24I) . i.e., as ?y* \ 0, is the gradient flow of 
Ent[o t ] on the Wasserstein space P 2 QR); see [9) ITOj. [42]. 

To conclude this subsection, we give an estimate on the second derivative y2,* ( t ) inside the bulk. 

Lemma A. 6 . Under the assumptions of Lemma \A.l\ the following holds. For ui £ [0, A], such that 
72,* (0) £ [E * — E/4, E* + E/4] we have 

|7Z*WI <CN S (1 + |72* ft)\) , 0 <t<C'N- 25 . (A.27) 


Proof. Let 0 < t < C'N 25 . For notational simplicity, we abbreviate here y w y = y2,* ft) and gt = Pt*■ 
We first compute 

^ [(Tft) (7«i,*)] = (Tft) (7 w,t) + [(Tp t )]' (7 w,t) iw,t 

= -i [T ((2 (Tet) (•) + -)Qt (•))]' (7 w ,t) - \ [T(T ft )]' (7 „,*) + M' ( 7 »,t) 7 , 

where we used (IA.24I) and (IA.14I) . (Here T((2 (Tp t ) (•) + -)g t (-))(x) denotes the Hilbert transform of 
the function y —>■ ((T g t ) (y) y)gt(y) evaluated at x.) Next, we note the identities 

y(T ft ) 2 -i^ = T ((T ft )ft) , 
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which follows from (IA.23I) and (T(g t ( •) •)) (x) = 1 + x (Tgt) (x), x £ M, which can be checked by hand. 
We hence obtain 


^ (Tgt) (7 w ,t) = ~ [(Tet) 2 ]' ( 7lM ) + \ [£]' ( 7w , t ) (T g t ) (.)]' ( lw , t ) 

+ [Tgt] (jw,t) iw,t + ■?f[£i] / (7w,t) j 

where we also used that (T(Tg t )) = — g t . Simplifying further and using (I A. 191) . we eventually obtain 


dt 


[(Tp t ) (7i«,t)] = [qV\ iiw,t) + ^[(Tp t )](7^,t) + 2[(Tet)]'('yw,t)'yw,t 


i r v ( \ 

+ ~2\£t\ W w,t) 2 " 


(Tgt)' (Tgt) 

Qt 


hw,t ) ■ 


We further compute 


d_ 

df 


(Tgt) 

gt 




d t (Tgt) 


gt 


(lw,t) - 


(Tgt) gt 




[(Tgt)]' 


gt 


7 t 


( T &) <?( 


ft 2 


(lw,t)'iw,t ■ (A.28) 


We next recall that we have the bounds |gt( 7 «j,t)| +1 (Tgt) ( 7 ™,t)| < C, as follows from the boundedness 
of mt (see Lemma ED, and 

|<9adm m t (x + ig*)| + I^Rem^x + i? 7 *)| < CN 5 , (A.29) 

for any x £ [E* — S/2, E* + S/2], see (IA.4I) . Thus, recalling that (Tg)/* (E) = Rem t (l? + i 77 *), 
7 rg/* ( E ) = Imr 7 it(£^ + i??*), we find 

I (Tgt)' (7»,t)| + |gt(7i«,t)l < CN S , |ft(7u,,t)| < CN S . 

Hence, differentiating (I A. 191) with respect to t and using (IA.28I) . (IA.29I) . we can bound 


\lw,t\ < ^l7w,t| + C ( 1 + 77* + 


V* 


V* 


{l + CN s + N s \^ Wtt \) . 


et(lw,t ) (gt(7m,t)) 2 , 

Finally, using that gt(lw,o) > c/2, for 0 < t < CN ~ 2S , as follows from the estimates (I A. 161) and (IA.17I) . 
and our assumption g( 7 uj,o) > c > 0, we immediately (IA.27I) . □ 


Corollary A.7. Under the assumptions of Lemma \A.6\ the following holds true. Let w,wo £ [0,A] 
such that 7ii)(0),7u, o (0) £ [E* — S/4,S* + S/4]. Then we have the estimates 

\ill(t)\<C, ItS'oWI < CN S , (A.30) 

and 

\i 2*(*) - ill (*)l < + c v *, |7S* 0 it) - ill (0)1 < cnH , (A.si) 

uniformly in 0 < t < CN~ 2S , with constants depending only on S, g and E and S. 

Proof. Since 7 kj( 0 ), 7 „, o ( 0 ) S [A* — S/4, A* + S/4], we have by Corollary IA.2I that j w {t), 7 „, 0 (t) ££ 
[A* — S/2, A* + S/2], for t < CN~ 25 , in particular we have gt(jw(t)), gt(7tu 0 W) > c > 0 for such t. 

Recalling the identity (Tg/ n *) (A) = R enit 0 (E + ig*) (as follows from (|A. 13|) and (12.51) 1. and the 
estimates on nit, d z mt derived in Lemma [A. II we conclude from that from (I A. 191) and (1A.20I) that 

l7t«o(*)l < c , \il*{t) -i2;(*)l < CN~ 1+6 \w - w 0 ] + CV , 

for all 0 < t < CN~ 2S . We further get from (IA.27I) that 

172; Ml < cn s , |72; (t) 72; (o) I < cnH . 

This proves (1A.31I) . □ 
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A.1.2. Proof of Lemma \4-2\ and of Lemma \4.3\ 

Proof of Lemma \4-2\ Without lost of generality, we can assume that ti = 0. Fix N £ N. From 
Lemma [A] TJ we directly get | gt(x) — p(f)| <C CN s t. Thus for t < C'N 26 , we get the first estimate 
in (14.81) . Next, we note that first and second derivative of m((z), z = E + i? 7 , are bounded on T>s /2 
by (IA.4I) . Thus the first derivative converges uniformly as 77 \ 0, E £ [A 1 * —E/2, A* + E/2], and we have 
Trd E Qt{E) = lim^o lmmt(E+bj). In particular, we obtain from (IA.4I) the second estimate in (14.81) . □ 


Proof of Lemma \4-3\ Without lost of generality, we can assume that t\ = 0 here. Let 77 * > 0 as 
in (IA.11I) . We then recall that we have | g?* (x) — Qt\ < for all x £ [A* — E, A* + E] and 

0 < t < C'N~ 2S . This follows directly from the definition of the Poisson kernel in (12.21) and Lemma lA.il 
Thus, using the definition of 7 in (IA.15I) . we must have I 7 /* (t) — 7 * (t) | < CN s r] t , i £ I a . By 
Assumption (1) of Theorem 12.11 mo( 2 ) extends to a continuous function on T>y,- Thus reasoning as 
in the proof of Lemma IA.ll we conclude that 7714 ( 2 ) extends to a continuous function on T>y ./2 for 
0 < t < C'N~ 2S . Hence, considering now 77 * as a free parameter (not depending on N) and taking 
77 * \ 0, we conclude from (IA. 191) that 


lim 

r 7 .N 0 




Qt(y)dy 

y-7i(t) 


7i(t) 

2 


i £ Ia ■, 


(A.32) 


for all i £ I G and 0 < t < C'N~ 25 . Here, we also used that > 0, gt(. 7 i(t )) > 0. Further, 

since 7 /*(f) converges uniformly to 7 j(t) for all 0 < t < C'N~ 2S , we can also exchange derivative and 
limit on the left side of (|A.32|) and we obtain (|4.11|). In particular, we have lim^ NOTTO) = 7 i(t), 
i £ I a , 0 < t < CN~ 2S . Thus (l4T3ll follows from (TOTll . 

Now we show (14.121) . For any fixed t we define the “continuous” quantiles 


r7«(t) 


Qt{x)dx = — 


u£ [0,1V], 


(A.33) 


and also the “half-quantiles” 7 j(i) := 7 i_i/ 2 (t) for any integer i. Since t is fixed throughout the proof, 
we drop the t argument. From (IA.33I) we get the regularity of the continuous quantiles: 


dy„ 

du 


1 

Ngi'yu) 


0(N ~ X ), 


in the bulk regime, where we used (14.81) . 
Setting j = £(i) for brevity, we can write 


dV g'( 7u) o(N ~ 2+s 1 

du 2 N 2 g{ lu y K ’ ’ 


f gt(y)d y 

pj-aN pj+aN + 1 

/ + + 

JrV- 7m 

J—OO Jjj-aN Jjj+frN + l 


q{v) dy 

V ~~ Tj 


(A.34) 


(A.35) 


The first integral can be written as (with 70 = — 00 and using (IA.33I) ! 

+1 ^ ^ ! | 3 ~44 P k+1 ( 7 fc - y)g{y)dy + P 1 g{y) dy 

to 4 V- Ti N ti Tfc -Ti ti 4 {V ~ 7j){lk ~ 7j) J-ooV-lj' 

(A.36) 


The last term is 0(N x ). The error term in the middle is bounded by 

J ~y 4 [ %+1 (t*- v)g(y)dy e { y )d y 1 

h 4 (y- Ti)(Tfc-Ti) -N ^ (Tfe-T,) 2 " ^ ti ( k-j) 2 ~ 


The third integral in (IA.35I) is estimated similarly. Finally, for the second integral we write 


n+° N+ ' g{y)dy y ( l 

h-s y~p k:l £l\=i V^Tk-Ti) 


P k+1 (Tfc - y) 0 (y)dy \ 
% (y - Ti)(Tfc - Tj )) 


P j+1 e{y)dy 

% y - Tj 


(A.37) 
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In the last integral we Taylor expand g(y) = g{lj) + 0(N 6 \^j — y\) by (14. 811 to get 

H i+1 g(y)dy P i+1 dy 


hi y - n 


eiv) i 

J'Y 


% y - 7j 


+ 0(N~ L ). 


Computing the integral explicitly and using 7 j = lj- 1 / 2 , 7 j+i = 7 j+i /2 we have 


rii+i 

hi 


d y 


y-io 


= log 


1 + 


7j— 1/2 - 2 7 j + 7j+i/2 


7j— 1/2 - lj 


= 0(N~ 1+S ). 


Here we used lj-i /2 — 7 j > c/lV and (IA.34I) to estimate the second order discrete derivative. 
Finally, we estimate the integral in the middle term in (IA.37II and we will show that 


aN 


E 

k:\k-j\=l ^ k 


7fe+1 ( 7 *: - y)e{y)dy 


(:V - lj)(lk - 7 j) 


= 0(N~ 1+S log N). 


(A.38) 


Clearly, (1A.38I) together with (1A.35I) . (1A.36I) . (1A.37I) and (1A.39I) imply (14.1211 . 
In the rest of the proof we show (|A.38I) . We write 


aN 

E 


f 7fe+1 ( 7 k - y)e(y)dy 


k:\k-j\=l ^ k 


(:y - lj)(lk - 7 3 ) 


(A.39) 


aN 

= E 

m =1 \ J 7j— 


(7 -j-m - y)e(y)dy 


0 y - io)(n-m - lj) 


L 


7j+m+1 (7j+m - l/)g(l/)dy 

7i +m (2/-7j)(7l+m -7j) y 


In the first integral, we replace g(y) with g(ij). From Taylor expansion, | g(y)~ g{lj)\ < CmN 1+5 , 


the error in this replacement is bounded by 

CmN 6 + 1 | 7 j_ m — y|dy CtoAT 5 p—+ 1 


IV 


'73-r, 


|y-7j||7i-m -7j| 


< 


dy 


IV 2 


< 


CN S 


’lj-n 


\v ~ IjWlj-m ~ lj\ Nm 


since |y j- m — lj I ~ m/N. We get a similar error when replacing g(y) with g(ij) in the second 
integral in (IA.39II . These errors, even after summation over m, are still of order 0(N~ 1+S logIV), 
hence negligible. Thus we get that (IA.39I) equals 


aN 

E 


Bin) 


f7j-m+l 


Ij-r 


'73- 


_ 1.j-m lj 

Next, using (IA.34I) . we have 

1 _ g( lo )N 


y-ij 


-d y 


eiij) 


lj+n 


lj 


Hi+m+l 


73+n 


7 j+n 


y 


y-n 


d y 1 + 0(N~ 1+6 log AT). 


0(N S ), 


lj+n 


Ij-m - lj rn 

so we get, after a change of variables, that (IA.39I) equals 


lj = 


TO 


g(lj)N 


+ 0{N~ 2+s m 2 ), 




rlj-m+i 


lj-m - y 


d y 


aN 


= -E 


N 


%-m V-lj 

flj-m + 1-1j-m 


/•73 + m + l ^ 

%+m V “ 7j , 

/*7j + m+l — 7j + 


7i+m ^dy ) + 0(A r_1+l5 log AT) 
u du 


m=l 


/Ti-m-73-m 7j “ 7j-m - U ’'73 + m—73 + rf 

The limits of integrations can be approximated as follows: 

1 


7i+m - 7j + u 


r yj—m r (j—m — 


2 Ngfrj) 


+ 0(mN~ 2+s ), 


lj+m Ij+m — 


2Ng( lj ) 


+ 0(N~ 1+S log AT). 


+ 0(mN ~ 2+s ), 
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1 


r )fj—m-\-1 'Jj—m — 


+ 0{mN~ 2+5 ), 


2N e ('y j ) ' " lj+m+1 lj+m ~ 2Ng{ lj ) 

Replacing these limits with their common values yields negligible errors, for example: 


1 


+ 0{mN ~ 2+s ). 


crN 

E 

m= 1 


_/\r nj-m+i-jj- 




u d u 


crN 


It i - jj-r 


- r<^E 

— 7 / ^ 


N mN s 1 1 


m N 2 N m/N 


= CN~ 1+S log AT 


Thus, with the notation d := l/2A r p(7j); we get 


tAT 


EMI = - E 


A 


udu 


udu 


+ 0(N~ 1+S log ./V). 


(A.40) 


—d 7 j Ij-vn u J —d 'Jj+m "fj + U 

Using that 7 ^ — lj- m = 7 j+ m — 7 j + 0(?n 2 A _2+l5 ) from (IA.34I) . we can write 

---=---+ 0{N 5 ), 

7 j+m - 7j- + w 7j - 7i-m + U 

for any u in the second integration regime, since here |u| < ^Nglfrl an< ^ 7 i+™ — 7 j T Ngfy-) + 
0(in 2 N~ 2+s ). Thus, replacing — 7 ^ with 7 ^ — 7 j_ m in the denominator of the second integral 

in (|A.40I) yields an error of order J] (A! /m)N~ 2+s = CN 1+S \ogN. After this replacement the two 
integrals in (IA.40I) cancel out: 


1 du 


I—d 7 j Ij — m u 


+ 


u du 


/ 


u du 


lj-lj-m+U J-dilj ~U‘ 


= 0 . 


We have shown that (IA.39I) is of order 0(N 1+l5 log AT), which finishes the proof of (IA.38I) . 


□ 


A.2. Remarks on Assumption (1) of Theorem 12.11 We conclude this Appendix with some re¬ 
marks on Assumption (1) of Theorem 12. II We consider the semicircular flow gt = J~t{o) started from g, 
for t > 0. As remarked earlier, the semicircle law g sc is invariant under the flow. It is then easy to 
check that m sc , the Stieltjes transform of g sc satisfies the bound in (12.151) for all t with 5 = 1. 

For many matrix models the distribution g, and hence also gt, are not explicit and checking As¬ 
sumption (1) directly may be not an easy task. In many situations, one can however use the smoothing 
effect of the semicircle flow to establish these estimates. The following example may be of some interest. 

Denote by C°’“(R), C' 0 ,a (C + ) the spaces of uniformly a-Holder continuous functions on R, C + . 
Assume that g £ C 0 , “(R), for some a > 0. Then the Stieltjes transform, mo, of g is in C 0 ,a (C + ). 
Adapting the proof of Lemma TA. II one can establish the following result. Abbreviate a t = l — e“L 

Lemma A. 8 . Assume that g £ C 0 ,Q (R). Then, mt{z), the Stieltjes transform of gt = Ft{(f), 
uniformly bounded for all z £ C + UR and t> 0. Moreover, there is a constant C, depending only on g, 
such that 


| m t {z) — nio(z)\ < Caf , 0 < t < 1, (A.41) 

and all z £ C + U R. Further, for all n £ N, there is C n such that we have the bounds 

\d?m t (z)\ < C n (a t lmm t (z )) a - n , t > 0, (A.42) 


for all z £ C + U R. 

Thus, running the semicircular flow from time t = 0 to time t\ = N~ Tl , t\ > 0, we see that 
Lemma [A. 81 implies the Assumption (1) of Theorem 12.11 for energies inside the “bulk” for the choice 
6 > (1 — a)ri. For the Wigner-like matrices of [2/5] typical choices for a are 1/3 or 1/2. 
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B. Persistent trailing of the DBM 


In this section, we prove that the time-dependent quantiles persistently trail the DBM up to a 

time-independent shift in the indices. More precisely, we have the following 

Proposition B.l. Consider a time interval [t 1; t 2 \ of length t 2 — t i = 0(N~ e ) with some small e > 25, 
where 5, given in (12.151) . is the regularity exponent of the initial data of the quantiles. Let A (t) be the 
solution of ( 12 . 111 ) and let 7 (f) be given by (12.91) . Suppose that 

p{|A i (t)-A j (t)|< iV€| ^~ j| , i , je ! a }< N ~ D , (B.l) 


for any D. Fix an index L in the bulk and let C(L) such that 

|A L(ti)--y m (ti)\<CN- 1+ t. 


(B.2) 


Then in the probability space of the Brownian motions {_B,(f) : i £ , t £ [ti, ^ 2 ]} we have 

f( sup |A L (f)- T ^)(i)| <CN~ 1+2 t) > 1-N-t. (B.3) 

\te[ti,t 2 ] / 

Notice that 7 e(L)(t) is a deterministic trajectory. This result therefore also shows that the typical 
fluctuation of the solution of the DBM is much smaller than the white noise term in (12.111) naively 
indicates. Indeed, the variance of the integral of this term is 


E 



t 2 — ti 
N 


which would indicate a behavior |A L(t 2 ) — Az,(fi)| > ( t 2 — tiY^N 1//2 . This is much larger than the 
actual value |A_Lp 2 ) — Ai(ti)| < CN ~ 1+ P 


Proof. Let 


Vi{t) := A i(t) - 7 ^ (i )(t) ■ 

Subtracting the DBM (12.111) from (14.121) and localizing it for the indices i £ I, we get 


dvi(t) 



d Bi(t) - ^2 - v j)dt - Wm d t + Ki(t ) d t , 

jei 


i £ I , 


with (time-dependent) coefficients 


Ba = 


1 

N\ A~ 


1 

a »)( 7 j - 7 %y 


ij e I, 


>v.=pi: 

k<jti 


1 1 

N (Afc - Aj)( 7 fc 


7 iY 


i £ I , 


and a deterministic error term |«j(i)| < N 

By (IB. II) and the spacing of the quantiles in the bulk, we know that 




b 


Wi(a) > 


b 

\j-L\- K\ + l ’ 


with b '■= N 1 * uniformly in time s £ [t -\, t 2 ] with very high probability. 
Let (s, t) be the propagator of the parabolic equation 


d Uj(t) 

d t 


y, - uj) - WiUi, 

jei 


i€ I , 


(B.4) 


(B.5) 
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then 


Vi(t)=Vi(t i)+ f y ^\fy{s,t)]j 
Jtl A 


/3 N 


d Bj(s) + Kj(s)ds 


(B.6) 


Since the propagator is a contraction in the supremum norm, the n(s) error term, after integration, 
gives a negligible error at most C(t 2 — t±)N~ 1+s < CN _1 . To estimate the main term, notice that the 
propagator depends on the sigma algebra £ s := {{Bi(u)}i^i : u € (s,f]} and is independent of the 
white noise at time s. Therefore 

EM*) - Vi(t i)| 2 = -J^E J* Y, + 0(N~ 2 ) . (B.7) 


Fix i £ I, s and t and define Wj := [f</ (s,t)]ij , which is the same as [& by symmetry. Then for 

any v > 0, we have 


< imcmhioo < 


CN* 




ll w l|i+^ 


by the heat kernel estimate on the Equation (IB.5I) : see Proposition 9.4 in [ 25 ] (the conditions of this 
proposition are guaranteed by (IB.4[) 1. By the L p -contraction of the semigroup for any p > 1, we have 


|w||i +v = \\^/(s,t)5i\\ 1+l/ < ||5i||i +J , = 1. 


Thus we get 


EM*) 


M*i)| 2 < 


CN$ r* 
N 1+ TT7 J tl (i 


^-^ds + 0{N~ 2 ) < C v 

S) ! + - 


_/y?+2y 

N 2 


< C 


N 2 ^ 

N 2 


after choosing v = £/2. Using Doob martingale inequality, one also has 


(B.8) 


E sup | Vi(t) - Vi(h)\ 2 < CE\vi(t 2 ) - Vi(ti)\ 2 < C-— . 

t<t 2 N- 

Setting i = L and using Markov inequality and combining it with assumption (IB.21) . we get (IB.31) . □ 
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